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Abstract 


The present thesis provides algorithms for some problems in computational geometry 
using prune and search technique. We solve the intersection radius problems in the 
plane for sets of line segments, rays, wedges, half planes and convex polygonal 
disks. Then we compute a centerpoint of a planar set of points. This is followed by 
description of a generalised technique of prune and search. 

The intersection raditts of a set of n geometrical objects is the raxlius of the smallest 
closed ball that intersects all the objects of the set. We have designed algorithms to 
optimally solve intersection radius problems for various kinds of objects. We show 
how the prune and search technique, coupled with the strategy of replacing a ray by 
a point or a line can be used to solve, in linear time, the intersection radius problem 
for a finite set of line segments in the plane. 

Next, the scope of this technique is enlarged by showing that it can also be used 
to find the intersection radius of a set of convex polygons in linear time. Moreover, 
it is immaterial if the set also contains other types of geometric objects like points, 
lines, rays, line segments, half planes and wedges. In fact, it is shown how such a 
mixed set of objects can be handled in a unified way; and this is the other important 
contribution of the thesis. Previously there existed no known algorithms to solve 
these intersection radius problems efficiently. 

The center of a set V of n points in plane is defined as the set of points c such 
that any closed half plane containing c contains at least [n/3] points of V. A 
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centerpoint is any point in the center. It can be viewed as a generalisation of the 
median of a set of reals. In the thesis, it is shown how a centerpoint of a finite 
set of points in the plane is computed in linear time. The described algorithm is 
optimal which significantly improves the O(nlog^n) complexity of the previously 
best known algorithm. We need to use suitable modifications of the ham-sandwich 
cut algorithm and the prune and search technique to achieve this improvement. 

The optimal sorting network by Ajtai et al. has a noteable feature that it approxi- 
mately sorts the input data in its intermediate steps. A technique is presented for 
solving computational geometry problems by exploiting this fact. This technique 
is a sequel to the Megiddo’s technique of designing serial algorithms by applying 
parallel computation algorithms. We get optimal linear time algorithms for some 
problems by applying this technique. We do this by the synthesis of prune and 
search and parametric searching. 
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Once there lived a village of creatures along the bottom of a great crystal 
river. Each creature in its own manner clung tightly to the twigs and rocks of 
the river bottom, for clinging was their way of life, and resisting the current 
what each had learned from birth. But one creature said at last, “I trust that 
the current knows where it is going. I shall let go, and let it take me where it 
will. Clinging, I shall die of boredom.” 

The other creatures laughed and said, “Fool! Let go, and that current you 
worship will throw you tumbled and smashed across the rocks, and you will 
die quicker than boredom!” 

But the one heeded them not, and taking a breath did let go, and at once was 
tumbled and smashed by the current across the rocks. Yet, in time, as the 
creature refused to cling again, the current lifted him free from the bottom, 
and he was bruised and hurt no more. 

And the creatures downstream, to whom he was a stranger, cried, “See a 
miracle! A creature like ourselves, yet he hies! See the Messiah, come to save 
us all!” And the one carried in the current said, “I am no more Messiah than 
you. The river delight to lift us free, if only we dare let go. Our true work is 
this voyage, fiis adventure.” 

But they cried the more, “Saviour!” all the while clinging to the rocks, making 
5 ? 

legends of a Saviour. 



Preface 


Geometry is one of the most ancient branches of mathematics. One of the earliest 
attempt to systematise Qeometry was made as early as in 300 B.C. by Euclid in 
Greece. This field also has distinction of being the first ever subject to profess 
of being rational and logical. At least it was believed by the people at that time 
and even now of standing on strictly formal grounds of logic. The whole system 
in Elements by Euclid is built on five basic geometrical postulates and numerous 
theorems are proved using only these few postulates. It can be said without any 
doubt that Geometry is one of the oldest field and scores of brilliant mathematician.*! 
and researches have worked in this field enriching geometry considerably by their 
discoveries and findings. 

Computational Geometry on the other hand is only of recent origin. It deads with the 
computational aspects, its theory and applications. The subject matter of Compu- 
tational Geometry ranges from theoretical computer science to that of more applied 
nature in Algorithms. To categorise the field might be difficult but it can be safely 
presumed that Computational Geometry is an offshoot of Applied Mathematics. It 
profusely uses the results of Combinatorial Geometry and constructs of Euclidean 
Geometry in its theoretical framework and the results of Analytical Geometry in 
its applications. When it deals with proofs of existences, theorems concerning 
Geometry, classifications of geometrical objects in different equivalent classes or 
enumerations of these then we are concerned with theoretical Computational Ge- 
ometry and when it deals with computing geometrical objects, or their attributes 
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or their count then we are concerned with application of Computational Geometry. 

It is important to know the types of problems in Computional Geometry to under- 
stand it better. The fundamental problems are that of construction of convex hulls, 
location of points in subdivisions, construction of Voronoi diagrams, etc. For a more 
complete list the pioneering doctoral work of Shamos can be referred. 

In this thesis we are concerned with the algorithmic aspects of Computational 
Geometry. We give algorithms to compute geometrical concepts such as intersection 
radius and center points. Intersection radius hcis got reference in the first ever 
collection of problems in Computational Geometry by Shamos where as the concept 
of center is dealt with in the book on geometry by Yaglom and Boltyanskii. The 
last part of the thesis deals with a generalised computational tool that can be used 
for geometrical optimisation problems. 

I wish to acknowledge my thesis supervisor Prof. Ashish Mukhopadhyay for his 
guidance during my Ph. D. work. I would like to acknowledge the work of researchers 
(especially Herbert Edelsbrunner and Nimrod Megiddo) before me in the field of 
computational geometry for providing motivation to persue this work. I thank my 
parents for leaving me alone to persue higher studies. I thank all my friends in and 
out of I.I.T. Kanpur for the company. I thank the I.I.T. Kanpur administration for 
bearing me. I thank the 
chance to do social work. 


workers’ cooperative of I.I.T. Kjunpur for allowing me a 
I finally thank myself for completing this thesis. 
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Chapter 1 


Introduction 


This thesis deals with a method of computing intersection radii for various kinds of 
geometrical objects and also of computing a centrepoint for a planar set of points. 
It also describes a general technique for applying the prune and search paradigm for 
solving geometric optimisation problems. 


1.1 Intersection Radius 


Suppose there arc several sites in a two-dimensional metric space and we have to 
choose a point p such that it is “nearest” to all the sites in question. Here the 
word “nearest” implies that the point p minimises the maximum distance of the 
sites from p. If a circle is drawn to “span” (cover) all the sites then it is called a 
spanning circle. When we draw a spanning circle centered on the “nearest” point p 
with radius equal to its distance from the furthest site(s), then this circle is called 
the minimum spanning circle. 

Phrases such as “spanning circle” and “minimum spanning circle” are meaningful so 
long as we are interested in spanning only point sites. These phrases are inadequate 
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to describe the problem of computing a smallest close<i disk that inters<‘fts 
objects such as lines, rays, line segments, etc. Thus the notion of stabbing is 
introduced. ^ 

When one object intersects or touches another it is said that the former stabs the 
latter. An object is called a stabber of a given set of objects if it stabs each object 
of this set. For example, the closure of the interior of the minimum spanning circU* 
above is the smallest closed disk that stabs the given set of points. 

Let C be a finite collection of geometrical objects in a <i-<iimensional hucliilean space. 
The stabbing problem consists of finding an object (the stabber) wliich inter.st'cts <'ach 
member of C. Typically, the stabber could be a line, a liyperplane or a disk, and C 
could be a collection of points, lines, line segments, rays, liyperspheres, po!yh<»drt)n 
or any mix of these. A survey of some recent results is available in tin* j>aper by 
Houle et al. [15]. 

A classical problem in this area is that of finding the intersection radius of a finite 
set of points in the plane, which is also known as the !-reniir problem (.'JO, .’Jl, (>j. 
It was shown by Megiddo, and subsequently by Dyer (22, 12, 11), how this can be 
solved in linear time for fixed d. However, until a recent paper {.‘Jj no attmnpt was 
made to extend this to a more complex set of objects than point.s or to a <-o!Ie<'tion 
containing different kinds of objects. 

Attempts have been also made recently to find more complicated .slabbers for th<!' 
stabbing problem, or to find the “best” stabber which optimises sonu* measurt* 
defined on the class of stabbers in question. Goodrich and Snoeyink [14] pre?sented 
an 0(n log n) algorithm to find a convex polygon whose boundary int<*rsects each 
of n parallel line segments. Rappaport and Meijer [24] shower! that a p<‘rin«der 
minimising polygonal disk that intersects each of n parallel line .segments can be 
found in 0(n log n) time. They have also extended their result to a set of isotiu'tic 
line segments. Mukhopadhyay and Kumar [27] have shown that for a set of parallel 
line segments an area minimising polygonal disk can also be found in O(rjlogn) 
time. Bhattacharya and Toussaint[4] gave an 0{n log^ n) algorithm for computing 
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the shortest line segment that intersects a set of n given line segments in the plane. 
This bound was later improved to 0(n log n) by Bhattacharya et al.[2]. 

In the present thesis, the stabber is a disk. The stabbing problem in this case is 
known as the intersection radius problem[13]. The intersection raditis of a planar 
collection of objects is the radius of the minimum stabbing disk of this collection. 
There is no known algorithm which solves the intersection radius problems for line 
segments. There are some variants of this problem which have been the subject 
of research of various researchers. As for example, when the stabber is a vertical 
line segment and the objects to be intersected are lines, the intersection radius 
problem is simply the Chcbyshev approximation of points in the dual space[5, 29]. 
This problem can easily be solved in linear time by transforming it into a linear 
programming problem. 

In this thesis, the intersection radius of a collection of line segments is computed 
by combining the prune- and- search strategy of Megiddo [21] with the strategy of 
replacing line segments with points or lines [3]. The scope of this technique is 
enlarged by showing that the intersection radius of a collection of convex polygons 
can also be computed in linear time. Further, it is immaterial if the collection 
contains a mix of other geometric objects such as lines, points, rays etc. for we show 
that it is possible to treat such a mixed collection of objects in a unified way. 


1.2 Centrepoint of a Planar Set of Points. 


We all have an intuitive idea as to what phrases like “the very center of the square” 
or “the very center of the city” mean. To capture this intuition in a quantitative 
way, the center of a set of n points, V, in is defined as the maximal subset of 
such that any closed halfspace intersecting this subset contains at least \n/{d -f 1)] 
points of V [35]. This subset is non-empty for any finite configuration of points (see, 
for example, [13]). Furthermore, it is closed and convex. A centrepoint is a member 
of the center of V. 



4 


On the real line a centrepoint is no other than a ine<lian of P . rinis a reiitrepoint 
can be viewed as a generalisation of the median of a set of reals. On the other haiui, 
the center can also be viewed as a particular fc-hnll of V. The k-huUo{ “P is a maximal 
subset (closed and convex) of 38*^ such that any closed halfspace inter-secling this 
subset contains at least k points of V. For instance, the 1-hull of V is its convex hull 
and the center is its fn/{rf -f l)')-hull. The property of halance<l partitioning makes 
the centrepoint useful for efficient divide and conquer algorithms in geometrical 
computing and large scale scientific computing[26, 32, 25, 36). Recently Donoho and 
Gasko have suggested that centrepoint can be used as “robust” and high “breakdown 
point” estimators for multivariate datasets [10]. 

The interesting algorithmic problem of computing a centrepoint has been considered 
by various researchers. Cole et al. gave an 0(n log® n) algorithm for computing a 
centrepoint of a planar set of points [9]. Subsequently, Cole improved this bound to 
O(nlog^n), using the powerful technique of slowing down a sorting network (8). In 
this paper, we propose an optimal linear time algorithm for computing a centrepoint 
of a planar set of points by using suitable modifications of the ham-sandwich cut 
algorithm for a pair of separable point sets [23] and the prune and search te< hni<pi<‘ 
of Megiddo [21]. 

Linear time algorithms, however, were known for computing an approritnatr or e- 
centrepoint [19, 32, 23]. We obtain this weaker type of centrepoint if wo decrease 
the lower bound, in the above definition of the center, to frjfl — e)/{d 4- 1 )], where 
0 < e < 1. Actually, Megiddo [23] only gave an algorithm for computing a partition 
of a (planar) set of n points with two lines such that each closed quadrant contains at 
least [n/4j points. An algorithm for computing an e-centrepoint, where 0 < e < 1/4, 
is implicit in this. 

The thesis proposes an optimal algorithm for computing a centrepoint of a planar 
set of points by using an interesting modification of Megiddo’s prtme-aii<i-.search 
technique[21]. This consists of adding a few extra points in each pruning ste[) .so that 
a subspace of the original solution space is retained, while ensuring a net deletion 
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of points. In the description of our algorithm, we a.ssume the usual RAM model of 
computation; the point set "P, however, is not assumed to be in general position. 


1.3 Prune and search in Slowed Down Sorting 
Networks 


A substantial part of Computational Geometry deals with designing efficient algo- 
rithms for given problems. More often then not, the solutions of these depend on 
the ingenuity of the inventor of the algorithm. However, there are also attempts 
to provide general tools to solve the problems of a specific type. Paradigms like 
divide and conquer, parametric searching, line sweep, prune and search, linear 
programming etc. go a substantial step towards this goal. 

The present thesis provides such a tool. We modify the method of parametric 
searching and slowed down sorting networks by Cole [8] to get a new technique. 

The parametric searching of Megiddo [20] is as follows — Let there be an efficient 
parallel algorithm for a problem A such that solution of A can be used in the solution 
of another problem B. Then, in some cases, we get an efficient sequential algorithm 
for B by exploiting the efficient parallel mechanism of the parallel algorithm for 
A. This technique has been applied to a wide variety of problems yielding efficient 
algorithms. In particular, we achieve good results for the parametrised problems 
that use parallel sorting algorithms in their solutions. We do this by replacing 
the evaluation of the parallel comparisons of an iteration in the parallel version 
by simultaneous resolution of these in the serial version. The running time of these 
algorithms is further improved by the introduction of weights in the comparisons [8]. 
We then simultaneously evaluate at least a fraction of total weight of these in every 
iteration to design more efficient algorithms. 
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We further improve the running time of the above algorithms by introducing, wher- 
ever applicable, prune and search in these. We do this by seeking to compute 
the fc-th largest element of the input set instead of seeking to sort it as in the 
previous techniques. For this, we run the sorting algorithm on a given input for a 
few iterations and then prune the set. It can be easily seen that this approach is 
useful only where pruning is applicable. 


1.4 Organisation of Thesis 


The organisation of thesis is as follows. In chapter 2 we solve the intersection radius 
problem for planar line segments in linear time. This algorithm is applied to all 
kinds of objects such as wedges, rays, halfplanes and convex polygonal disks in the 
next chapter. Moreover, it is shown how such a mixed set can be treated uniformly. 
In chapter 4, we present a linear time algorithm for computing a centrepoint of a 
finite planar set of points. We present a generalised tool for applying prune and 
search to a given problem in chapter 5. We conclude and discuss clirections for 
further research in the last chapter. 



Chapter 2 


Intersection Radius of a Set of 
Line Segments 


The intersection radius of a finite set of geometrical objects in a (^-dimensional 
Euclidean space, is the radius of the smallest closed ball that intersects all 
the objects of the set. In this chapter, we describe optimal algorithms for some 
intersection radius problems in the plane. We first briefly review a linear time 
algorithm by Megiddo to determine the smallest disk that intersects a finite set 
of points. Next wc show how the prune and search technique, coupled with the 
strategy of replacing a ray by a point or a line can be used to solve, in linear time, 
the intersection radius problem for a finite set of line segments in the plane. 


Previously, no efficient algorithm was known to solve the intersection radius problem 
for line segments in the plane. 
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Figure I: Spanning circles and stabbing disks of points and line segments 

2.1 Introduction 


Suppose there are several sites in a two-dimensional metric space and we have to 
choose a point p such that it is “nearest” to all the sites in question. Here the word 
“nearest” is to be understood in the following way: the nearest point p mininiise's 
the maximum distance from the sites to p. This problem frequently arises in practice 
in the installation of centralised facilities e.g. such as installation of a transmission 
centre. If a circle is drawn centred at p with radius equal to the distance from p to 
the farthest site, to span all the sites, then it is called a minimum spannmg circle 
of the sites (Fig. 1(a)). 

The phrase “spanning circle” is meaningful only if we are interested in spanning or 
covering points and is inadequate to describe the problem of computing the smallest 
closed disk that intersects other objects such as lines, rays, line segments, etc. 'rhus 
we introduce the notion of stabbing for these objects. 

When one object intersects or touches another we say that the former stabs the 
later. A stabber of a set of objects is an object that stabs each member of this set. 
For example, the closure of the interior of the minimum spanning circle above is the 
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smallest stabbing disk of the given set of points. 

'I'hc concepts of spanning circles and stabbing disks, though related, are different 
in several ways. To span a set we only need to consider extreme points of the 
individual members of the set. For example, to span a set of line segments we only 
need to compute a spanning circle for the endpoints of these. However, this does 
not help in the computation of minimum stabbing disk of line segments which is 
frequently smaller. Similarly, in general also, it does not help to characterise the 
minimum stabbing disk if the minimum spanning circle of the given set is computed 

(Fig- 1(b)). 

The intersection radius of a planar set of objects is the radius of the minimum 
stabbing disk of this set. Now this question naturally arises — how to compute the 
minimum stabbing disk/intersection radius for a given set of objects efficiently in 
linear time. This problem is progressively solved in this chapter and the next, first 
for simple objects, viz points and lines, and then for more complex ones. 

The problem of computing the intersection radius usually reduces to a non-linear 
programming problem except when the stabbed set consists of straight lines alone. 
Not only for these but also for their counterparts in higher dimensions, i.e. hyper- 
planes, the intersection radius problem similarly reduces to a linear programming 
problem [3]. For fixed d, this latter problem can be solved in linear time, using the 
algorithm of Dyer or Megiddo or Clarkson [11, 21, 7]. Or instead, a recent simpler 
and more efficient randomised algorithm of Seidel [28] can be used that runs in 
0{n • (d -h 1)1) expected time where n is the number of constraints in the linear 
programming problem. 

The case of points also has been fully tackled [21]. This was done by exploiting the 
fact that the bisectors of points are straight lines. These bisectors can be seen as 
linear functions and therefore the methods of linear programming problem had been 
applied to this case too. 


The organisation of this chapter is as follows. In Section 2.2 we give a few definitions 



10 


and discuss the motivation of the method to compute the intersection radius of points 
in linear time. The solution of the intersection radius problem for points in plane 
is presented in Section 2.3. In Section 2.4 we design algorithms to compute the 
intersection radius of lines and rays which are used to compute intersection radius 
of line segments in plane. 


2.2 Preliminaries 


Let P be a set of n points in the Euclidean plane aiul let r be the centre of a 
spemning circle of radius r. If d(c,p) denotes the distance of c from a point p in P, 
then clearly, 

r > max d(c,p) 

Denoting by, (/(c,?^), the distance of c from the set V, which is also the ra<lius of 
the smallest spanning circle centred at c, we have 

d{c,V) = tmad(c,p) 

The intersection raditis, 1R{V), of V is the radius of the minimum spanning circle. 
Thus, 

1R{'P) = nund(c, P). 

We can similarly define the intersection radius of a set of lines, rays, line segments 
etc. We will denote minimum spanning circles by min-circle and minimum stabbing 
disks by min-disk . The centres of these in the unconstrained and constrained 
versions of the problem will be denoted by unconstrained centre and constrained 
centre respectively or just centre when the context is clear. 

Let us first see how intersection radius of V is computed. Consider the min-circle 
of three points in the plane. The circumscribed circle of the triangle formed by 
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these points is also a spanning circle of these. However, it is not necessary that 
this circle is a minimum too, unless the triangle formed by the points is an acute or 
right-angled triangle. In the case of obtuse-angled triangle the min-circle is obtained 
by drawing a circle with its diameter as the largest side of the triangle. A naive 
algorithm of complexity O(n^) to compute the intersection radius of V now is same 
as computing the radius of the largest of min-circles of all subsets of V of size three. 
It can be shown that the largest of these min-circles spans V. 

A better algorithm, that runs in (9(n^) time, can be designed as an improvement to 
this naive algorithm by fixing two points in each iteration such that these two points 
are on the circumference of the spanning circles. These points are then dynamically 
updated until the min-circle of V is finally obtained. The details are as follows. We 
start with two fixed points, such that the points in V lie on the same side 

of pq. The min-circles of three points are drawn for p, q and every other point of 
"P. Let the point corresponding to the largest of these circles, which incidently also 
spans V, be r. If i\pqr is obtuse-angled then we delete the obtuse-angled vertex 
and repeat the same steps with the largest side as chord, otherwise the minimum 
spanning circle of Apq'r is the min-circle of ^(Fig 2). 

We might think of improving this algorithm further by fixing only one point and 
then computing the spanning circles. But, since we need two more points to fix a 
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spanning circle, and we have a quadratic luunher of choices, some other approach 
is needed. We therefore consider fixing a line on which the centres of the spanning 
circles are constrained to lie. Further, we also have to think of some way of not 
considering all of the quadratic choices that determine a spanning circle. This leads 
us to the important idea of prune and search. 

The radius of the min-circle of T such that centres are ronstraiiu'd to he any |)oint 
p is a convex function of p, f{p). In particttlar, when the domain of this function 
is restricted to a line, /(•) is still convex. Therefore hy examiinng the gradient of 
/{•) at any arbitrary point p on a cpiery line we can tell whetln'r f[p) is miiiiinutii 
and if /(p) is not minimum then on wliich side of p on th<* iin<‘ the ininitnum of 
/{•) is attained. We will see in the next section how these observations h<*lp us iu 
computing intersection radius of points in linear time. 


2.3 Intersection Radius of Points in the Plane 


When the objects to l)e intersected an* only points, tin* algoritlmi of Dyer or Megiddo 
[12, 22] for the l-centre problem can be used. A brief account of tiie algorithm is 
presented in this section, such that, it can be further gtnierali.sed to include other 
kinds of objects. 

Let V be the set of points whose int<>rs<*ction raditis is to b<* <!elermiin*<! and l<*t its 
size be Np. 

To begin with the following question is answered: Which are tin* points that can h<* 
pruned. To answer this, we note that min-circle is determined by a set of only two 
or three points. These points are farthest from the centre of min-circh* among all 
the points in T. We try to determine these by removing the other points of P in an 
organised manner as follows: 


IJ the centre of min-circle is contained in a half plane determined by the 
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Figure 3: Localisation of Centre of Min-disk of V 

bisector of some two points of V then the point which is neai'cr to the 
centre can he dropped from V in the computation of intersection radius. 


To prune away a constant fraction of points according to this rule, we need to localise 
the centre of min-circle with respect to a constant fraction of bisectors. The crux of 
the method is — How can this be done? 

We further make note of the fact that, given a non optimum spanning circle of 
radius r, the radius of min-circle is smaller than r and therefore, the centre of min- 
circle will be at a distance smaller than r from any of the points in P. This has 
the following consequences. First, for any point P, we can localise the centre in 
a small conical region whose apex is P. This can be done by first computing the 
smallest possible spanning circle centred at P and then computing intersection of 
interior of circles with radii equal to the radius of this circle centred at points on 
its circumference, of which at most three need to be considered for the purpose at 
hand (Fig. 3). Secondly, for any line L, we can localise the centre in one of the sides 
of L, by similarly computing the min-circle whose centre lies on L first and then 
computing the conical region, which can be shown not to intersect L, with respect 
to the centre of this circle. 

To compute intersection radius of P in linear time we consider the following related 
subproblems. Subproblem 2 and Subproblem 3 can be solved using oracles of 



Siibprohlcin 1 aiui Sul)i>rol>l(‘ni 2 rosix't livcly. 'I’hc a)ijv<*X!ty of fuiirlioii 

is crucial in both the solutions as it also provides the direction, consequently a half 
space, in which the search for the centre of inin-circle is localised. 

Subproblem 1 Compute the min-circle ofV such that the centres of the spanning 
chclcs are at a fixed point /\ 

Subproblem 2 Compute the min-rircle of 'P such that the n nires of ih< spanning 
circles lie on a fixed line L. 

Subproblem 3 Compute the jnin-ciir.lc of 'P. 

The solution to the first subproblem is the simplest. We compute the maximum of 
the distances of points in V to the fixed point P and choose this distanc<* as the 
radius of min-circlc. This computation takes linear tim<'. 

The min-circle, whose centre is constrained to lie on line L, is determine<i as follows. 
We pair the points arbitrarily in [A^p/2J pairs and ctunpute the p«‘r{M*ndirular 
bisectors of each pair. If any bisector is parallel to th<‘ line L then wo drop tin* 
point that is nearer to L in the corresponding pair. In case the bisector is same as 
line L then any one of the points in the pair is dropjx'd. After this, Subprtdjlem I 
is solved for the median point, P, of the intersection points of L with tlu* rest of 
bisectors. If P is the constrained centre then the computation is over, otherwise 
one point corresponding to each of these bisectors can be dropped by localising the 
constrained centre on L with respect to P. Thus it can be seen that at least j 
points are dropped in one iteration. We iterate until fewer than four points remain 
when any straight forward algorithm is applied. This is the techuKpje of prune and 
search in which the size of input set is truncated by a constant fraction in each 
iteration. It is easily seen that it runs in linear time in the present case. 

The solution of Subproblem 3 is obtained as follows. We pair the points in P and 
form their perpendicular bisectors. I^^et the median slope of non-vertical bisectors of 
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these be Sm. The bisectors of slope strictly larger than .s,„ are [jain’cl with bisectors 
of slope strictly smaller than s^- As it can be seen, this pairing may not cover all 
the bisectors, if the number of positive slope bisectors is not equal to the number 
of negative slope ones. Next, intersection points of each of these bisector pairs are 
computed. 


First we localise the centre of min-circle with respect to a vertical line which divides 
the total of the vertical bisectors and the intersection points above in two halves. If 
the centre of min-circle lies on this vertical line then the inin-circle is determined, 
otherwise let Jy be the halfplane, determined by this line, in which the centre of 
'min-circle is localised. Next we localise the centre of min-circle with respect to a 
line with slope such that it divides the total set of the intersection points in Jy 
and the bisectors with slope in the above. As previously, if the centre of min- 
circle lies on this line then we terminate our computation, otherwise let Jm be the 
halfplane determined by this line in which the centre of min-circle is localised. Then 
wo drop a point which is nearer to the centre of min-circle for each corresponding 
bisector that does not intersect JyC\Jm. ■ 


It can be easily shown that there are at least L^p/4j bisectors which are one 
of these — the paired bisectors, the bisectors with median slope or the vertical 
bisectors. Thus, in every iteration at least [A^p/16j points axe dropped. These steps 
are repeated with the truncated set of points until there are fewer than sixteen 
points when any straight forward algorithm to compute min-circle is applied. This 
algorithm runs in linear time. 


As we will show later in the next section, the method of computation of intersection 
radius of a set of lint's, £, also uses the same principles. 
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2.4 Intersection Radius Problem for Lines, Rays 
and Line Segments in the Plane 


In the previous section we presented a linear time algorithm to determine min-circle 
of a set of points V. In this section we shall design an algorithm to d<*t<*rniiine the 
minimum stabbing disk, min-difik , which intersects a given set of line segments, 
rays and lines. 

The real difficulty to apply the same technique to compute the intiu sectiou ra<!ius of 
line segments/rays is that the bi.scctors of these can lx* curved bin's. nature of 
these depend on the relative positions of the line segments/rays in the |>lan<*. The.se 
may contain parts of bisector of two endpoints (a straight line), parts of bisector 
of an endpoint and a supporting line (a parabola) and a part of bise<-lor of the 
supporting lines of both (an orthogonal pair of straight lines). It is the presence 
of parabolas in the set that does not allow us to apply the pruning techni(ju<* to 
these bisectors. However, the curved parts, namely parabolas, in the bi.sectors of 
line segments owe their existence to the bisectors of Iin<* and f)oints only. So, we 
can alternatively look at a similar problem in ord<T to circumvent this dilliculty. 
This problem is — How can we determine the inters<-ction railiu.s of a mixed .set 
of objects in the plane containing both lines and points. An aildixi int<‘rest in this 
latter problem is that its solution is also used in the solution of th<' former prtdilem. 

Further, each line segment is considered to be eipiivalent to two oppositely directecl 
rays having the endpoints of the line segments as their respective tails. 'I'hus the 
intersection radius problem for a set of A//., line segiiKuits can !><‘ r<*<!iieed to the 
intersection radius problem for a set of 2Nig rays in the plane. I'he latter problem is 
solved by combining the prune and search strategy of Megiddo [21] with the novel 
idea of replacing a ray by a line or a point in the pruning st<‘p. 'I’liis is tione !)y 
breaking up the problem into the following two subprobleins: 


Subproblem 4 Given a set of lines and points in the plane, compute the smallest 
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7'adius disk that intersects these. 

Subgroblem 5 Given a set of rays, points and lines in ike plane, show how a 
fraction of the rays can be replaced by lines or points such that the intersection 
radius of the new set is the same as that of the original one. 

We will use the subscripts l,p and r to denote functions related to lines, points and 
rays, symbols L, P and R to denote a line, a point and a ray and denote finite sets of 
these by C,V and TZ, respectively. The functions gi{x,y),gp{x,y) and gr{x,y) have 
the following definitions; 

5/(-T,2/) = rnax{di(/v,(x,j/)) I /, € £}; 
gp{x,y) = max{dp(P,(x,y)) P € p}; 
gr{x,y) = max{dr(P,(x,y)) ReH}, 

where di(L, {x, y)), dp{P, (x, j/)) and dr'{R, {x,y)) respectively denote the (Euclidean) 
distance of the point (x,y) from a line L, a point P and a ray R. We define g{x,y) 
as 

g{x,y) = max{ 5 j(x,j/), gp{x,y), 

It can be easily shown that the functions < 7 ;(x, y),gp{x, y) and gr{x, y) are all convex 
so that g(x,y) is also convex. 

Let S be any finite set of lines, points and rays. The convexity of this last function 
enables us to answer two key questions about the minimum radius disk that inter- 
sects S. The first is that if we constrain the centre of the minimum radius stabbing 
disk to lie on a fixed line L then to which side of a given point (a,^) on L does the 
constrained centre lie. The second is that given a line L to which side of it does the 
centre of the minimum radius stabbing disk lie. Without any loss of generality, we 
can take L to be the x-axis of an orthogonal frame of reference. 
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Lot UR answer the first question. Clearly, wo can coinpiito ,9(0,0) for tlie wt 5 in 
0 {Ns) time, where Ns is the size of S. Let S' be the subset of objects who.so 
distance from (a, 0 ) is <?(«, 0 ) or, more simply, those ti»at touch the disk. Since 
9(0,0) is convex, if the contact points of all the objects in S' lie to the !eft( right) of 
the vertical line through (o,0) then the centre of the constrained minimum radius 
stabbing disk lies to the left{right) of (o, 0 ). Otherwise, (o, 0 ) itself is the required 
centre. 

The second question is also easily answered. We compute a ininiinum ra<liu.s stabbing 
disk whose centre is constrained to lie on the line /,. If fix' «'(>i)ta< t points of (he 
ol)j<x;ts in «S' .span an arc greater than or e<jual to a setni ciK h* of the disk buiindary, 
then the computed disk is the required minimum raditis spanning disk. Otherwise, 
the centre lies in the same halfspace, <let.ennined by L, as in wlii< b tlx- mid point of 
the chord of the above spanning arc lies. This again follows from the convexity of 
the function g{x^y). 

We first present a linear time algorithm to determine a minimum radius stabbing 
disk that intersects a given set of points and lines. When the obj<‘cls to be inter.sect<xi 
are only points, the algorithm of Dyer or Megiddo (H, 21] for the I centre problem 
can be used. The situation is more complex when lin<!« are also iuclixh'tl it! the set of 
objects to be intersected. Our algorithm uses the basic prune and search techni<jue 
of Dyer or Megiddo, referred to above. We describe the techni<|ue for Uiu*8 only. W<* 
will then argue at the end of .section that the ad<lition of point.s does not change tlx‘ 
underlying concept. 

2.4.1 Intersection Radius Problem for Points and Lines in 
the Plane 


Here we first present a linear time algorithm to <let(‘rinin(‘ tniii-disk of a set of bin's, 
C. Let the size of C be Ni. The min-disk of C can be computed by transforming 
the problem to a linear programming problem and then solving it using any of the 
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known methods [3]. The algorithm given hei-e uses a different approach and applies 
the basic prune and search technique of Dyer and Megiddo [11, 21]. 

As in previous section, we progressively solve three similar problems. These problems 
are — First, computation of the stabbing disk of minimum radius when the centres 
of stabbing disks arc at a fixed point P, computation of minimum radius stabbing 
disk when the centre of stabbing disks lie on a fixed line L and lastly, computation 
of the unconstrained niin-disk. 

The first of these problems is solved in linear time by computing the maximum of 
the distances from each line in £ to point P. The required min-disk has this distance 
as its radius. 

In the solution of the other two, as in the case of points, the corresponding bisectors 
of pairs of lines play an important role. These bisectors are either orthogonal pair 
of angular bisectors (for pair of non parallel lines) or parallel midway line (for pair 
of parallel lines). The lines are dropped in the computation of intersection radius 
by the application of th<^ following rule: 


If the centre of min-disk is localised in one of the regions, defined by the 
bisector(s) of two lines in the set C, then the line that is nearer to the 
localised region can be dropped. 


To prune away a constant fraction of lines in C by applying this rule we need to 
localise the centre of the min-disk with respect to the bisector(s) of a constant 
fraction of pairs of lines. The localisation of centre with respect to a single line is 
done by computing the function g{x, y) as mentioned above. We discuss next how the 
localisation of the centre of min-disk, with respect to the bisectors (angular/parallel 
bisectors) of a constant fraction of pairs of lines in £, in linear time, is done. 
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I A CouHtrnhiod Version of (ho Prohivui 

We discuss the problem of computing the min-disk for a set of Hut's, £, siu h that 
the centre of this disk is constrained to lie on a line L, which is assum<*d here to he 
the x-axis of an orthogonal frame of reference. We also assume that £ ha.s at most 
two lines parallel to L, since any other lying hetwwu tlu'se two <ioes not change the 
intersection radius. Furthermore, out of these two lim*s we can prune away the one 
that is closer to L. So there is no loss of gen<*ra!ity in a.ssuining that then* i^ at most 
one line parallel to L. To the remaining /V/(> {Nt — \ )) !ii»*s we apply the prtme 
and search strategy in the following way. 

We first identify a group of roughly half the line-pairs such that one of their hi.seclors 
does not intersect an interval on L which contains the ce'ntre of the constraim'd min- 
disk. This can be done in linear time as follows. 

Wo compute the n»odian ;r,„ of th<' intersections of all the lines with I'lien we 
determine the value of at this point and use this to locate the conslrainetl 

centre on L with respect to (smiO). If the constrained centre, <!enoted by {x'.O), 
lies to the right (left) of this median point, the required half consists of tho.Ht> Hues 
whose intersections lie to the left (right) of the median. We label the Hih's /,,, with 
I b- will suffice to discuss the case in which tin* coa.slrained cj*ntre 

lies to the right. 

Consider the set of line-pairs (/y2,_i, £2,), with 1 < ? < 

For each parallel line-pair it is clear that we can {>run<' away th<* on<' that is chwer 
to Xrn- The pruning mechanism is non-trivial for non-parallel pains, since we have 
to do a finer location of the constrained centre. Of these, we can prune away the 
line closer to Xm if one of the angular bisectors is parallel to line £. 

Each non-parallel pair (£21-11 £21) an associated pair cjf angh' hiseefots. 'I’he 
intersection with L of one of these does not lie between tin; interse<‘tions of t h<* lin<s< 
themselves (Fig. 4 ). Let d; be this intersection point. 



21 



Figure 4: d,- does not lie between the intersections of lines L^i and L 2 i-i- 
dm lies to the left of x*. 



As before, we locate the constrained centre with respect to the median dm of those 
intersections among these that lie to the right of x„i- Either x* — dm, or we have 
the following cases: 


Case 1: x* < dm 

Consider a d, that lies to the right of dm- Since x* lies to the right of all 
the Li's, the angle bisector that intersects L at d,-, its associated pair, 
together with the lines that these bisect give rise to the configuration 
of Fig. 5. Since x* lies as shown we can prune away one of the lines of 
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the pair, (/> 2 ,-i, /-a,), it can ho easily shown that we can prune roughly 
one-eighth of the liiu’s we started with. 

Case 2: x* > dm 

It is easy to see that in this case also approximately one-eighth of the 
lines we started with can be thrown away (Fig. 4). | 


The above pruning takes 0{Ni) time. We repeat this process until no in<>r<’ liiuvs can 
be pruned, when the optimal .solution ran be obtained l»y some brute for< «' met hod. 
It can be shown that no more than eight liiu's are left at thi.s stag*'. If tin* initial .s<‘t 
of lines contained a line parallel to L, the radius of the constrained min-di.sk is the 
maximum of the optimum value obtained and the distance from this parallel line to 
L. The total running time of this algorithm is easily se<'n to be 0[Nt). 

If S also contains points, each pruning step is carric<l out in two substc-ps. In the 
first substep we prune points, followed by lines in the .second or the other way rouml. 
To prune points, we first pair them arbitrarily, and compute the median point of the 
intersections of L with the bisectors of these pairs. 'Phen we rh'tertnine the min disk, 
£>, centred at the median point, for the set S. Next, to determine on whi<'h si<ie 
of the centr<? of I) on b the coii.straiiu'd <<'ntrc Hi's, wv <'xaiiiin<' how tin* poiuls of 
tangencies of lines of S, that touch the disk D, and the points of that lie mi its 
circumference, are distributed with respect to the vertical line through the <'entr<‘ of 
D. When we prune lines next, again points of S are ignored similarly in the stejjs 
in which we determine Xm and dm- 

In summary, we note that when pruning objects of one kind, the objects of the other 
kind become transparent whenever we need to find a point on i to serve as the c<*ntre 
of stabbing disks for all the objects in tlie set currently umler consideration. 

If Ni and Np be the number of lines and points in S respectively, it is easy to .see that 
we prune away at least one-eighth of the total number of object.s, viz Nt 4 Np. We 
repeat this process until we cannot prune any more obj<.*cts. The constraiiu'd centre 
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Figure 6: One of the bisectors of Lj and Lj [B-j in tlio figure) does not intersect the interior of 
LL. 


can then be determined by a brute force algorithm. The whole process requires 
linear time. 


I The Unconstrained Centre Problem 

Assume that £ is the set of lines in 5. We compute the min-disk of £ in the 
unconstrained case. A.s before, we will indicate later how to handle the addition of 
points. 

We pair up the lines in £ arbitrarily and compute the angle bisectors of each pair. 
In the degenerate case of a pair of lines being parallel, the angle bisectors reduce to 
a single line parallel to and equidistant from the lines that make up the pair. When 
a pair of lines have distinct angle bisectors, these make up an associated pair. 

The following observations are crucial. Given a pair of intersecting lines, if we can 
locate the region containing the centre in a quadrant defined by the angle bisectors 
of the pair then we can prune away one of the lines. For a pair of parallel lines we 
can do the same if we can locate this region in a halfplane determined by their “zero 
angle bisector”. We indicate below how we do this for a fraction of such pairs. 
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(’onsidor first th<‘ subset of vertical bisectors. We locale the unrou8traine<l centre 
with respect to a median vertical bist'ctor Li- Assume that the centre of tnliedisk 
lies in the left halfspace of this median line. Clearly about half of the vertical 
bisectors do not intersect Jj. We now compute the median slope of the non- vertical 
bisectors and pair up arbitrarily in this subset a bisector which has slope greater 
than the median slope with one which has a smaller slope. The bisectors in each of 
these pairs necessarily intersect as they have unequal slopes. Next we compute the 
median of all the unpaire<l bisectors whose slopes are e(jual to the nx'dian slope. We 
again locate the unconstrained centre with respect to this median bisector I. j, whicli 
has median slope. Assume that the centre lies in the halfplane ./^ below tins !im>. 
Clearly about half of the median slope bisectors that lie above tlo not intersect 
J 2 - The centre, now lies in J\ H 'I'o be able to prune any lines, we ntvd to refine 
the location of tin.' centre still birtlu'r as in the worst ca.se, w<* may n<*ilher have any 
vertical bisectors nor any whose slopes are equal to the median slope. 'I'lu'refore, we 
do the following with respect to the above pairs of intersecting bisector.s. 

We first locate the unconstrained centre with respect to a vertical line through the 
median of the x-coordinates of their intersections. We can a.ssume, without any lo«.s 
of generality, that the unconstrained centre lies to the left of this lijie. We now 
project, parallel to the median slope, the intersection point.s that lie t<» right of thi.s 
line onto the y-axis. We then locate the unconstraiiuxl c<‘ntre with respect to a line, 
parallel to the median slope, passing through the median of these projections. Again 
there would be no loss in generality if we assume that tin* unconstrained centre li<*8 
below this line so that now it lies in the lower left quadrant <ietermined by this 
line and the earlier one. 


This ensures that at least a fourth of the pairs of bisectors in this class have their 
intersections in the upper right quadrant V R. Consider one such pair. Hln<-<* the 
bisector with slope smaller than the median slope does not intersect LL, at least 
one-eighth of the bisectors whose slopes are not equal to the median slope <k> not 
intersect LL. We note that this argument does not depend upon which quadrant the 
unconstrained centre lies in. Thus, we have a set of bisectors which do not interswt 
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Figure 7: The region J = LLCi JiD J 2 . 


the region J = LL fl Ji fl J 2 in which the centre of min-disk lies (Fig 7). In the 
worst case even this set of bisectors does not enable us to prune any lines because 
there may be no associated pairs of bisectors or bisectors due to a pair of parallel 
lines in this set. Further, in the worst case, it may happen that all of the associated 
bisectors of this bunch intersect J. 

We need to do one more final refinement of the location of the unconstrained centre 
with respect to these associated bisectors that intersect J. If we repeat the above 
steps with the above associated set of bisectors we get a region J\ containing the 
centre and a fraction of these bisectors which do not intersect it. Each of these, 
together with its associated pair from the earlier set of bisectors we found that do 
not intersect J, contain in one of its quadrants the region J fl J'. Region J fl J' 
contains the unconstrained centre. Therefore one of the lines whose angle bisectors 
these are can be pruned. At the same time we can prune away one of the lines of 
each parallel pair whose “zero angle bisector” does not intersect J D J'. 

Calculation shows that at least [iV//64j of the lines are pruned away. 

We repeat this process until no more lines can be discarded. It can be shown that 
there are no more than 64 lines in £ at this stage and we use some brute force 
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Figure 8: Ray R is replaced by the supporting line of If is replaced by its tail point 
method to compute the unconstrainct} centre. 

When points are also included in the set, S, of objects to be inter.sertetl, we gt» 
about the pruning step in exactly the same way as in the con.strained case. In t‘ach 
such step we throw away a well determined fraction of the point.s and the lines. 
Repeating this process we get a linear time algorithm for the unconstrainerl centre 
problem for 5. 

2.4.2 Intersection radius for rays 

The constrained problem for a set TZ of Nr rays in the plane can be rediiceti to the 
problem of pruning for a set of rays, lines and points as detailerl below. 

For each ray Rj € R, consider the line //,■ normal to it and passing through its 
tail. We compute the median x„, of the inter.section.s of these normals w’ith the 
constraint line L and locate to which side of this nH‘dian on L the constrained 
centre lies. Suppose the constrained centre lies to the right of this median. Then 
for each normal which intersects L to the left of this median point we replace t he 
corresponding ray by a line or point according to the following replacement rule : 

If the ray and the median point lie in the same halfspaee of the (wo 
halfspaces determined by the normal then the ray is replaced by the line 
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which cojitams the ray. Othertoise, the ray is replaced by its tail point 
(Fig. 8). 


It is important to note that neither of these replacements changes the radius of the 
relevant stabbing disks of the original set of rays. Thus, if Nr be the number of rays 
in 71, at least [fVr/2j rays are replaced by either a point or a line and our new set 
of objects consists therefore of lines, points and rays. Next, all the points and lines 
are considered for pruning as described in the § 2.4.1. We thus discard a fraction of 
the rays from further consideration. These two substeps are iterated until no more 
objects can be disc.ard<'d, when any brute force method can be applied to compute 
the min-disk centred on L. It is easy to check that the algorithm runs in linear time. 

lo solve the unconstrained version of the intersection radius pi'oblem for a set of 
rays, we need to replace a fraction of the rays by points or lines in linear time. This 
can be done as follows. 

As in the constrained case, we start with the normals through the tails of the rays. 
Proceeding identically as in the case of the unconstrained problem for lines, we 
determine a region J which contains the unconstrained centre and is not intersected 
by at least one-eighth of these normals. We do not need to iterate twice as in the 
case of lines. Only one iteration suffices for the replacement of a fraction of rays. 
The ray corresponding to each of these normals can therefore be replaced by a line 
or a point according to a similar replacement rule. 


If the ray and localised region J lie in the same halfspace of the two 
halfspaces determined by the normal then the ray is replaced by the line 
which contains the ray. Otherwise, the ray is replaced by its tail point. 


Thus in a single iteration we replace about one-eighth of the rays by lines or points. 
Next We use the two-step pruning process on the set of lines and points generated so 
far to throw away a fraction of them. Repeating these two substeps on the modified 
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set of rays, lines and points, and then applying any forn* inetiiiMl whrii iiu 

pruning takes place we get a linear time algorithm h»r inlei'ic* tiuii ptohU-m of hue 
segments in the plane. 


2.5 Concluding Remarks 

In this chapter we have described optimal algorithms for rojupiiting the •mmUtrst 
radius disk which intersects a set of line segments in the plane using a novel apptoai h. 
It would be worth investigating whether a similar a|)pioa« h < aii l»e Uit other 
kinds of stabbing problems. 



Chapter 3 


An Optimal Algorithm for the 
Intersection Radius of a Set of 
Convex Polygons 


In the laist chapter it was shown how the minimum stabbing disks and intersec- 
tion radii for finite sets of points, lines, rays and line segments can be computed. 
These were computed in linear time by combining the prune and search strategy 
of Megiddo [21] with the strategy of replacing line segments and rays by lines or 
points [3]. In this chapter, we enlarge the scope of this technique by showing that 
it can also be used to compute the intersection radius of a finite set of convex 
polygons in linear time. Moreover, it is immaterial if the set also contains other 
types of geometric objects such as points, lines, rays, line segments, half planes and 
wedges. In fact, we will show how to handle such a mixed set of objects in a unified 
way. 

The bisectors of line segments, as we have seen in the previous chapter, are made 
up of lines and (parabolic) curves. Here too the bisectors of convex polygonal disks 
are formed by not only parabolic curves but also regions of non-zero width. Such a 
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case arises wh<‘ii a pair of rotivex polygonal disks have non-empty interseetion and 
we have these regions of non zero finite area included in the bisector(8) of theae. We 
present a method in this chapter in which this non linearity of bisectors does not pos«* 
any problem and we compute the minimum stabbing disk of a finite set of convex 
polygonal disks in linear time. But before doing this we first formalise the notions 
of replacement, substitution and localisation. 'I'hese concepts were implied when we 
used these in their specific contexts in the previous c!iai>ter hut lun-ertheless, were 
not discussed there. 


3.1 Introduction 


Let C be a finite collection of objects in the Euclidean plane. Tlie stabbing [irohleni 
consists of computing an object (the stabber) which intersects each member of C. 
Typically, the stabber could be a line, a disk, a line segment etc. ami C ctudd he a 
collection of points, lines, line segments, rays, circles, polygons or any mix of tlie.se. 
A survey of some recent results is available in the paper hy iloule et a!, jlh). 

The intersection radius of a finite colUrtion of geometrical ohj«*cts is the radiu.H of 
the smallest clo.sed ball that intersects all the objects in the collirlion. A « 
problem in this area is that of finding the intersection radius of a finite set of 
points in the plane, which is also known as the I-centr<* problem (30, 31, fij. 'flie 
corresponding disk is called the minimum stabbing <lisk. It was shown hy Megidilo, 
and subsequently by Dyer [22, 12], how this can be solved in linear tinw. More 
recently, Welzl [34] has given a randomised algorithm for this problem which runs 
in expected linear time. However, until an earlier paper hy Bhattacharya et a!. [3] 
no attempt was made to extend this to a more complex colle<'tiot» of objects than 
points or to a collection containing <lifl'<Tent kiiuls of ohj<'ct.s. 

New attempts have been made recently to find iiiori* complicated .stabh<*rs for 
the stabbing problem, or to find the best stabber which optimises some mea.sure 
defined on the class of stabbers in question. Goodrich and Siioeyiiik [14] pr<*sented 
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an 0(n log ri) algorithm to find a convex polygon whose boundary intersects each 
of n parallel line segments. Rappaport and Meijer [24] showed that a perimeter 
minimising convex polygon that intersects each of n parallel line segments can be 
found in 0(n log n) time. They have also extended their result to a set of isothetic 
line segments. Mukhopadhyay et al. [27] have shown that for a set of parallel line 
segments an area minimising convex polygonal disk can also be found in 0(n log n) 
time. Bhattacharya et al. [2] gave an 0{n logn) algorithm for computing the shortest 
line segment that intersects a set of n line segments in the plane. 

Bhattacharya et al.[3] showed that when C is a collection of line segments the 
intersection radius can be found by combining the prune and search strategy of 
Megiddo [21] with the strategy of replacing line segments with points or lines. This 
was discussed in the last chapter. In this chapter we enlarge the scope of this 
technique by showing that the intersection radius can also be found in linear time 
when C is a collection of convex polygons. Really, it is immaterial if C also contains 
other geometric objects like lines, points, rays etc. We show how it is possible to 
treat such a mixed collection of objects in a unified way. 

The organisation of this chapter is as follows. Section 3.2 contains the necessary 
geometric and algorithmic preliminaries. In Section 3.3 we describe the algorithm 
and analyse it in the following section. The last section contains conclusions and 
directions for further research. 


3.2 Preliminaries 


In the rest of this chapter we will adopt the following notation scheme: collections 
of objects will be denoted by letters in script style such as C,£,P, ..., objects by 
capital letters such as L,P,R ,.. ., and points by small letters p,q,r,... Let the size 
of a collection C be denoted by Nc . 

The algorithm that we shall describe in the following section is based on three 
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important terhniquos. 'I'lM'sr arc: n plact uti nl of a unuplcx incitniciiH al u!«)ci t liy 
a set of simpler objects, localisation of the rciitrc <jf juinimum jiiabiMUR iti a 
more restricted region and thereafter filtmtwn «>f some <if the ohjeiis whnh do not 
play a role in determining the centre of the nnniinnm stahhing disk logethi'r. the 
latter two make up the prune anti search techni<|ne that was iirsl fornialtHed in 
Megiddo [21]. 

We discuss these strategies, replacement, lorahsattan and filinitwn, in detad in the 
next three sections. 

Let C be a collection of n objects in the plant* where ea« li ofij*'*'* is eithei a |><iint, a 
line, a line segment, a ray, a wedge or a convex ptdygtnial tlisk. The jutddem is to 
determine the minimum stabbing tlisk, min tlisk, of the t nlbs lion C. 

In the following discussion J is a fixed subset of plant* containing the centre t>f iiiiti- 
disk. In the general case, J is initially the whole plane, and gets sumller as the 
computation proceeds. 

In the constrained case, J is always a subset of the line on width the .entreit of 
stabbing disks are constrained to lit*. All the tliscussion in this se. tion is also 
applicable to the constrained case whi*re J is such reslrictetl to a line. 

Let C and C denote an object and a collection of objects respectively. 'I he ilistance 

of a point p from <7, d(p, C), is the shortest Kuclitlean di.stance of p frtmi ( ’Hie 

distance of p from C, d(p,C), however, is the Iarge.st of all dislanecs from p to the 

objects in C. The stabbinff mdius, SR[J,C), where the centre of the slabbing tlisks 

are constrained to lie in J, is the minimum of the distances from }><jints in J to C, 

The intersection radius, IR(C), is the unconstrained stabbing radius. In snimiiary. 
wc have 


d{p,C) = infd(p, 9 ), 
d(p,C) = maxd(p,6-), 
SR(J',C) = min d(p,C), 


M.A 
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Figure 9: A polygonal disk, and its substitution by wedges (one of the wedges seen slightly shifted) 

IR(C) = SR(3^^C), 
where rf(p, 9 ) is the distance from p to q. 

The distance bisector of two objects S and T is defined as the set of all points for 
which the distance to S is equal to the distance to T, i.e., it is the set 

{ P I dip, S) = d(p, T), p € 3?^} . 

3.2.1 Object Replacement 

To simplify the computation it is helpful to divide a convex polygonal disk into a 
set of elementary parts such that the distance function is not modified. This allows 
us treat these in a uniform manner in the algorithm. To see how we can divide a 
convex polygonal disk, let us first see how the distance from a point p in to a 
•convex polygonal disk is computed. This distance is either equal to the distance 
from p to some vertex, or it is equal to the perpendicular distance to some side, 
or it is zero when the point lies inside the disk. A wedge is the non-reflex region 
bounded by its two infinite sides. It can be easily seen therefore that the distance 
from p to this disk is the maximum of the distances to the wedges formed by taking 
the vertices and their adjacent sides (Fig 9). So if a convex polygonal disk in C is 
substituted by these wedges then the intersection radius of the resulting collection 
does not change. 
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Figure 10: A line segment, and its Rubstitution hy twn rays (Rrrn 


Each line segment can he similarly replaretl hy two <>j»jiositely <!iriH't«*ii tays whiwe 
intersection is the line segment itself (Fig 10). 'I'he (lisfanee of tlii> Ime •H-ginetit fruin 
a point p in is cither equal to the dislanre of the endpoints or it is equal l«i the 
perpendicular distance to the supporting line; this clistanee i.s clearly equal iti the 
maximum of tiie distances of p from the two rays. 

These examples motivate the following theorem. 


Theorem 3.1 Let C be an object and V a collection of objtcfii. nuch that tUp, ( ') v 
for all p€ J. Then for any collection C. 

SU{./,CU{f’}) = SR{./.C’UT>), 


Proof: From the definitions of ijiUusection raditis am! di.statjce fuitctions, we have 

SR(J,CU{C}) = nnn max{d(j?,C),cf(p,(')} 

- nwn max{d(p.C),d(p,I»)} 

= min d{p,CUV) 

= SR(J,CUX>). 

This proves the theorem. | 


When C is a line segment, a ray, a convex polygonal <lisk, a halfjdane t»r a werlge, 
we can apply Theorem 3.1 to replace C hy some simpler ohject{8). 
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Ray is replaced ^ Ray is replaced 

by tail point ^ by supporting line 

\ 

Figure 11: A ray it, its associated normal and regions 

In the two examples mentioned above J is the whole plane. This allows unconditional 
replacement of line segments and convex polygonal disks by rays and wedges, respec- 
tively. Let us take another set of examples, when J is only a proper subset of the 
plane, C is a ray, a halfplane or a wedge, and the replacement is done conditionally. 
Let J? be a ray in C. The normal to this ray through its tail point divides the plane 
in two halfplanes. If we localise the centre of min-disk in one of these halfplanes 
then the distance from any point of the localised region to the ray is equal to either 
the distance to the tail point or the perpendicular distance to the supporting line 
of the ray R (Fig. 11). This means that we can replace R either by its tail point or 
its supporting line. 

Similarly if J is localised in the interior of a halfplane, H, in C then the distance 
from any point of J to if is zero, otherwise if J is localised in the interior of H then 
the distance from any point of J to if is equal to the distance to the boundary of 
if. We can respectively discard if or replace H by its boundary line in C if these 
cases arise (Fig 12). 

Let us draw an outward normal to each of the sides at the apex of a wedge W. 
The whole plane is then divided into four (unequal) quadrants by these normals and 
sides of W. If we somehow localise the region J in one of these quadrants then the 
distance of the points in J from the wedge has one of the following values: 0, when 
J is localised inside the wedge; perpendicular distance to one of the sides, when 




36 



Figure 12: A balfpl&ne H, its associated boundary line and rrgmna 


EepiacemeM 

by 


Deletion of I Reptacenwni 
Wmlfs / fay Hide 

/ Wed*e 

Figure 13: A wedge W, its sMociated normals and trgtons 


J is localised in the region bounded by a side and the norma! to it; or, distance 
to the apex, when J is localised in the region bounded by two normals. We can 
respectively discard W, replace W by the relevant side, or replace W by the apex, 
in these cases(Fig 13). 


We shall henceforth view a convex m-|on as the collection of m wedges defined by 
the vertices and the sides incident on them. Our problem, therefore, is equivalent to 
that of computing a min-disk for a set of N„{= £ m, sum taken over all the convex 
polygons) wedges. Likewise, we shall view a line segment aa a pair of oppositely 
directed rays defined by its end points [3]. 
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Figure 14: Stabbing disk is a min-disk 


3.2.2 Localisation of the Centre of Min-Disk 

In the discussion below we assume that the objects are convex sets of points so that 
all distance functions are convex. 

An object C is said to touch a disk of radius r centred at point p if r = d(p, C) , 
and the vectors from p to the points of contact are said to be its contact vectors. 
The following theorem characterises min-disk. 

Theorem 3.2 A stabbing disk of radius r centred at point p is min-disk iff all of 
its k contact vectors, r,-, with I < i < k, are linearly dependent satisfying 

53A,ri = 0, 

j=i 

where Xi > 0, for 1 < i < fc, with some A,- ^ 0 and k >2. 

Proof: The proof makes use of the behavior of the distance function in the 

neighbourhood of p, which depends only on the contact vectors of the stabbing disk. 

Let the contact vectors r,- satisfy the relationship given in the theorem. Then for 
any arbitrary vector v, we have 

Z; A<(ri • v) = 0 

i<k 
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hence, (r^ • v) < 0 for some i, with 1 < « < A*. ‘hn ohjrn i hr f ‘ NW. 
if p is displaced to by a small amount in the dirrrtjun of v and .f ^v' it ihr 
corresponding displacement of the new contact point. p[, id ihr ohjn t ( * ( I- ig H) 
then the square of the new distance, of the new centre /»' fr..in the ..hjr. i ( ‘ „ 

df = Ir; - (5v - ^v')P 

= jf.P + l^v - ^v'p - 2(r, ■ ^v) + 2(r, /*v') 

> Ir.f 

> r^. 

Since C is convex and r< is normal to it. r, • is non negatiw and this justifies 
the inequality above. Since d{p,C) is the maximum of all distances fmin p to the 
objects in the collection C, d!(.,C) will also increase in the diret imn v I hus 

d(p',C) > < 

> r 

Hence, as v is an arbitrary vector, p is a local minimuin of the function rf( .t") Since 
d(.,C) is convex, p is a global minimum too. 

For the proof of the converse, we assume that the said disk is the intn disk. Then wc 
prove that the centre p lies in the relative interior of the convex hull of the contact 
points. 

We will prove this by contradiction. We assume the contrary that the centre of min* 
disk lies does not lie inside the convex hull of contact points. Then we c«i compute 
a line which separates the point p and the convex hull. I.»et v be the normal vector 
to this line, contained in the same halfplane as the convex hull. I'hen v ■ r, > 0 for 

all t, with 1 < t < *. FVom this it can be seen that d(p,C) strictly decreases along v, 
since 
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Figure 15: Centre of Min-Disk lies inside convex hull of cont 2 ict points 

< 

< |r,p — • (2r, — ^v) 

< |r.f 

< r? 

< d{p,Cf. 


This means that the radius of the min-disk is not minimum. Since this cannot be, 
therefore the assumption that the centre of the min-disk lies outside the convex hull 
of contact points is incorrect. 

Let the vectors to p and the contact points be p and p,-, with 1 < * < fc, respectively. 
Now, an interior point of a convex hull can be written as a positive linear combination 
of the extreme points. Thus 

P = S 

i<k 


%<k 



40 


where A,- > 0, with 1 < t < fc. Hence, 


- P) = 0 

«<* 


and therefore 


53 A,r, = 0. 


I 


The above theorem is also applicable in the constrained case with a slight modifi- 
cation that here the contact vextors projected onto the constraint line are positive 
linearly dependent. Once we have computed a stabbing disk, if it is not a minimum 
one as determined by Theorem 3.2, we can further localise the region in which the 
centre of min-disk lies with the help of the theorem below. 


Theorem 3.3 Given a stabbing disk of C centred at p, with k contact vectors r, , 
with 1 < i < k, the centre of min-disk lies in the set J' given by 

where Hi is the halfspace, normai to r„ that passes through p and contains the 
contact point and the tangential object. Furthermore, 

SR(J,C) = SR(J',C). 


Proof: Let the radius of stabbing disk be r. We have to show that centre of 
min-disk lies in every Hi. Let the corresponding object to //, be C. The radius of 
min-disk has to be smaller than or equal to r. Therefore the centre of min-disk is 
at a distance smaller than or equal to r from every object in C and in particular 
from C. The feasible region of this centre will thus be a subset of i/, (because C 



41 



Figure 16: Localisation of Centre of Min-Disk 

is convex). Every point not in Hi is at a distance larger than r from the object C 
(Fig 16). This proves the theorem. | 


We can use this theorem to localise the centre of min-disk with respect to any 
arbitrary fixed line, L, in the plane. We first compute the minimum stabbing disk 
such that its centre is constrained to lie on L and then compute the region given 
by the expression of the theorem with this disk as the reference. Since the region 
contains only those points where the value of stabbing radius is smaller than the 
radius of given disk, therefore it necessarily does not contain any point of L. Thus 
we get a region that is fully contained in one of the halfplanes determined by L. 
Further, if we localise the centre of min-disk to a region J' we then onwards replace 
and filter objects with respect to newly located region J' and not the earlier located 
region J. 


3.2.3 Pruning or Filtering Objects in C 


We now need a suitable criterion to discard objects from the collection C that are 
irrelevant in determining the centre of min-disk. The following theorem provides 
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Figure 17: Filtering line* and point* 


such a criterion. 


Theorem 3.4 Lei there he a pair of objects in C such that region J is contained in 
one of the regions determined hy their distance hiseciorfsj. Then one of these objects, 
whichever is nearer to J, can be discarded from C without affecting the inff rjirrhon 
radius. 


Proof; FVom the definition of the distance bisector, if the region in J i» contained 
in one of the regions defined by the bisector then, for every point p 'm J, distancr of 
p from one of the objects is sdways smaller than the distance from the other object. 
Let tins nearer object be C. Since the intersection radius is the maxtinuin of the 
distances of objects of C over the points of the region J, the object C , its distance 
being smaller than the distance of the other object in the pair, does not play a 
role in the determination of min-disk. Hence, the object C can be deleted from the 
collection C. % 


For example, refer to Fig. 17 above. The angular bisectors of two lines ij and ia 
and the perpendicular bisector of two points Pj and P3 are shown. The nearer of the 
lines, Lif is deleted because region J (the shaded region in the figure) is contained 
in the one of the four quadrants defined by the angular bisectors of h and 1%. 


\ 
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Similarly, the nearer of the points, Pi, is deleted because the region J is contained 
in one of the halfplanes defined by the bisector of Pi and Pj. 


3.3 Intersection Radius Problem for Convex Poly- 
gons 


In this section we present a linear time algorithm for computing the min-disk which 
intersects a collection C of convex polygons. In fact, we will be more general and 
assume that C also contains points, lines, rays, line segments, halfplanes and wedges. 
Since each line segment can be considered to be equivalent to two oppositely directed 
rays having the end points of the line segments as their respective tails, and a convex 
polygon as a collection of as many wedges as vertices (refer § 3.2.1), we will assume 
that C consists of lines, points, rays, halfplanes and wedges only. 

We first solve a constrained version of the problem and then the unconstrained one. 
The motivation behind this is as follows. The intersection radius is the minimum 
of all stabbing radii over the domain of the plane. This problem can be formulated 
as a non-lineax programming problem. The usual method of solution of such a 
problem in consists of locating the optimal point with respect to a hyperplane 
in one dimension lower, by first solving the problem restricted to this hyperplane, 
and then locating the optimal point by computing the gradient of the minimising 
function. 

Furthermore, to replace or filter objects, we need to locate the centre of min-disk 
with respect to either bisectors of pairs of points, bisectors of pairs of lines, the 
wedges or the normals to rays. So we apply the same technique, as in non-linear 
programming, to localise the region containing the centre. We note that we need 
not do this for every object in C in a single step. It suffices to do this only for a 
fraction of these objects. So, we successively localise the centre of min-disk with 
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Figure 18: Localisation of Centre of Min Disk in Linear T ifiie 


respect to a group of lines, which axe chtMen in such a way that there is always a 
fraction of the objects in C which are either discarded or replactni. 

First we provide an efficient linear time implementation of the construction suggested 
by Theorem 3.3. Then we solve the constrained and unconstrained rases of the 
intersection radius problem in successive subsections. 


3.3.1 An Efficient Implementation of Theorem 3.3 

Suppose we axe given a stabbing disk centred at c. We have to locate the centre of 
min-disk in a region, J, such that c lies on its boundary. We do this in the following 
way. 

If the stabbing disk does not touch any object then its radius can be shrunk until 
it touches at least one object. This can be done in linear time by computing 
the maximum of the distance from c to the objects in C. Suppose there are k 
^contact vectors of the stabbing disk with k>l. If we compute the intersection of k 
halfplanes as in Theorem 3.3 by constructing the convex hull in the dual plane, then 
we will need 0{k logk) time. This can take 0{Ne log Nc) time in the worst case. In 
our application, however, this intersection can be computed in linear time because 
all these halfplanes have point c on their respective boundaries. Computation of 
their intersection is same as the convex hull computation in the dual plane of a 
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set of collinear points. We only need to compute the two extreme halfplanes and 
the optimal centre lies in the intersection of these two halfplanes which is a wedge 
(Fig 18). If the intersection of the halfplanes is a single point, namely c, then c 
is the centre of min-disk. In the constrained case we need to locate the centre of 
minimum stabbing disk on a line L. For this we compute the part of the line L that 
lies inside this wedge. 


3.3.2 The Constrained Centre Problem 


In the constrained version of the problem we determine a min-disk whose centre 
is constrained to lie on a line L. Let T, £, 72., H and W be disjoint subsets of 
C containing points, lines, rays, halfplanes and wedges respectively such that their 
union is C. We show how we can process points, lines, rays, haJfplanes and wedges 
separately. We then unify these into a single iteration of the algorithm. 


Points: We pair up the points in V arbitrarily and compute the median of the 
intersections of their bisectors with the line L. The centre is then localised 
with respect to this point on one of the half lines. We filter out at least one 
fourth of the points corresponding to the bisectors that do not intersect the 
localised region. We repeat this step until no filtering takes place. 

Lines: We pair up the lines in £ arbitrarily and compute their distance bisectors 
(angular bisectors for non-parallel lines, and mid-way line for parallel lines). 
We divide the intersections of these bisectors with the line L equally into four 
intervals on L. The centre of min-disk is localised in one of the intervals by 
doing a binary search on the boundary points. We filter at least one fourth 
of the lines corresponding to the distance bisectors that do not intersect the 
localised region. We repeat this step until no filtering takes place. 


With rays, halfplanes and wedges we proceed in a slightly different way. 



46 



Figure 19: Intersection of objecU with J m the cpi«lr«oe4 pH>hWm 

Rays: We draw normals at the tail points of the rays in K and nimpulr ihr »i»«lian 
of the intersections of these normals with the line i. The cenlrr of itmi disk is 
then localised with respect to this point. We can replace a ray. whi«w normal 
does not intersect the localised region, cither by its tail |K»mt or Us supporlmg 
line. We replace at least one half of the rays in each itrraiiun Wr rr|H-al this 
step until no replacement takes place. 

Halfplanes: We localise the centre of min-disk on the line L with res|»ect to the 
median of the intersection of line I with the boundaries of halfplanes in H. 
We replace a hadfplane in H by its boundary line or discard it if its iHUindary 
line does not intersect the localised region. We replace/filtrr at least one half 
of the halfplanes in cac^ iteration. We iterate until no filtration/ replacement 
is done. 

Wedges: For each wedge in W, we draw outward normals to the sides at its apex. 
We divide the inter8ectk>n(8) of the normals and the sides with the line L 
equally into four inter's^s of L. We then localise the centre of min-disk to 
one of these intervals. Ail those wedges (at least one fourth) whose sides and 
normals do not intersect the localised region sure either filtered out or replaced 
by a line or point. We iterate until there is no replacement or deletion. 

It is easy to design an 0{NclogNc) algorithm, based on the above facts. We 
first convert all the rays in C to lines or points, then we filter or replace all the 
halfplan^ and wedges in C, and finally compute the intersection radius of a collection 
containing only lines and points. To obtain a linear time algorithm, we have to treat 
the objects in a more unified way in each iteration of the algorithm. 
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The basic idea is as follows: The derived lines obtained as above from different 
objects of C are classified as either active or inactive, depending on the fact, whether 
region J which contains the centre of min-disk is localised with respect to these. 
Further, active lines among these are assigned following weights: derived lines of 
wedges with four and three intersections with J - 16; derived lines of wedges and 
line pairs with two intersections with J - 18; derived lines of rays and line-pairs with 
one intersection with J - 27; and derived lines of point-pairs with one intersection 
with J - 36. At the beginning of computation all the derived lines that intersect 
L are active. In each iteration we first localise the centre of min-disk on L with 
respect to the weighted median of intersections of line L with the active derived 
lines. The locating of the centre of min-disk is done by invoking Theorem 3.2 and 
Theorem 3.3 to determine on which side of the weighted median intersection point 
the constrained centre lies. Finally, we replace and filter those objects, none of whose 
derived lines are active, by invoking Theorem 3.1 and Theorem 3.4. We repeat this 
until no active derived lines axe made inactive. 

It can be easily seen that in every iteration half of the active lines of about half 
the total weight are inactivated. New derived lines of at most only three-fourth of 
this discarded weight are added. So there in net reduction of weight by a constant 
fraction. Hence the algorithm is linear in total weight. 

The details are given in the algorithm below. V is the active derived set of lines. 

Algorithm 1 CONSTRAINED-CENTRE 

Input: Bisectors of point -pairs and line-pairs, normails to rays, 
boundeury lines of half planes and sides of emd normals to wedges. 
Output: Centre of constrained min-disk 
begin 

T> < — input set 
do 

Compute the weighted median intersection, P , of X> with L 
Locate J with respect to P on L that contains the 
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co&8trai&«d castra 
Update t> 

Filter aad replace relevant objectt of the collection C 
Update £ 

while is updated 
enddo 

Detemine the constrained centre b; eoM brute force Method 
end I 


The proof of correctness of the algorithm CONSTRAIHED-CEHTRE i« a direct rtiiier- 
quence of the theorems in Section 3.2. 


3.3.3 The Unconst rained Centre Problem 

The tasks of filtration and replacement are more involved in the unronst rained case. 
We first describe a method by which we can localise the region that rontains the 
centre of min-disk which does not intersects a finite fraction of a given set of lines. 
As in the constrained case, we then show how we can process points, lines, rays, 

halfplanes and wedges separately. Finally, we unify these separate steps into a single 
iteration of the algorithm. 

Let LOCALISE be the procedure which, pven a set of lines £ with integral weights, 

determines a region J that contains the centre and is disjoint from at least one 

eighth of the total weight of lines. Also, our algorithm terminates if the centre of 

min-disk is found during a call to the procedure COSSTMIIfED-CEIITRE in LOCALISE. 

So we only need to consider the case in which the algorithm does not terminate in 
this manner. 

We first compute the weighted median slope, of the non*vertical lines in £ and 
hen divide the lines into three sets, £>, and £«, having slope greater, smaller 
than and equal to s,„ respectively. Let C, be the set of vertical lines in £. We pair 
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the lines in £> and £< arbitrarily using their weights. We then compute intersection 
points of these pairs. The weights of both members of the pair should be equal and 
this weight is assigned to the intersection point. For example, a line with weight 3 
can be paired twice — once with a line of weight 2 and second time with a line of 
weight 1 whereas the intersection points have weights 2 and 1 respectively. We can 
ensure the number of pairs to be linear as all the lines have integral weights only. 
Note that if the total weights of the sets £> and £< are unequal then some of the 
lines may remain unpaired. 

Next we compute a vertical line, Ly, which divides the total weight of the intersection 
points of the above pairs and the vertical lines in £y into two equal halves. We 
localise the centre of min-disk in one of the halfplanes determined by Ly, by using 
the solution of the constrained centre problem with Ly as the constraint line. Let 
Jy be the halfplane determined by Ly in which the centre of min-disk lies. 

We then compute a line Lm with slope Sm which divides the total weight of points 
and lines — intersections points of the above pairs which lie outside Jy, and lines of 
£= — into two equal halves. We again locate the constrained centre on this line and 
determine on which side of it the centre of min-disk lies. Let Jm be this hadfplane 
so that the centre of min-disk now lies in Jy fl Jm- 

It is easy to see that at least one eighth of the total weight of lines in £ are disjoint 
from Jy n Jto- These lines cross the opposite quadrant Jy D J^. 

The description of LOCALISE is as follows: 

Algorithm 2 LOCALISE 

Input: Set of Lines £ 

Output: Localised Region J 
begin 

Compute the intersections of an arbitrarily (equal weighing) 
paired set of lines in £, one of larger slope 
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than the weighted median slope and the other of tealler slope 

Localise Centre in Jy with respect to a vertical line which 
halves the total weight of intersection points and 
vertical lines in C 

Localise Centre in Jm with respect to a line of weighted Median 
slope which farther halves the total weight of intersection 
points in 7^ and weighted Median slope lines in C. 

J • Jfn n 

Return J 
end I 


Once we have localised the centre of min-disk, we repeatwUy apply Thmri’rn 3.1 
and Theorem 3.4 on wedges, halfplanes, rays, lines and points to prune/rrplarr the 
objects in the collection C. We discuss separately the cases of points, linm, rays, 
halfplanes and wedges. 

Let V, £, TJ, W and >V be the subsets of C containing all points, lines, rays and 
wedges in C respectively. ' 

Points. We pair up the points in P arbitrarily and compute the perpendicular 
bisectore of these pairs. We use procedure LOCALISE with this set of bisectors 
aud filter one point corresponding to every bisector that does not intersect 
the localised region. Per iteration this filters out at least i/,. of the ,»inls 

(Pig. 20(a): The centre is localised in Lower Left quadrant, /./. and the point, 
bisector crosses t/pper LHgbt quadrant, UR). 

Lines; Lines are filtered by localising the region containing the centre of min disk 
m a quadrant defined by the bisector(s) of a pair of lines. When the linm are 
paraUel the distance bisector is a line parallel to and equidi.t«,t from them. 
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We will call a pair of angular bisectors mates of each other in the following 
discussion. When we invoke procedure LOCALISE with this set of bisectors as 
input, there is no way to ensure that there exists at least a fraction of these 
bisector pairs which do not intersect LL. So, we have to do the localisation 
twice in the following way. 

We first pair the set of lines in C arbitrarily and compute the bisector(s) 
of these pairs. We localise the region containing the centre of min-disk with 
respect to this set of bisectors. Let LL be this localised region. We consider all 
the bisectors that do not intersect LL. In the worst case, all these are angular 
bisectors such that their mates intersect LL. We invoke LOCALISE again with 
these intersecting mates as input. A fraction of these mate-bisectors do not 
intersect the localised region, say LL', returned by this second invocation of 
LOCALISE. These together with there mates from the previous invocation of 
LOCALISE do not intersect region LL fl LL' which contains the centre of min- 
disk. We filter one of the lines of each pair of lines whose distance bisector(s) 
does not intersect LL H LL'. This filters out V 64 of the lines in C (Fig. 20(b)). 

Rays: Rays in % axe replaced by localising the region containing the centre of min- 
disk in one of the halfplanes defined by the normal at the tail point of a ray. 
We do this by invoking LOCALISE with the set of normals at tail points of rays 
as input. This replaces at least i/s of the rays in "R, (Fig. 20(c)). 

Halfplanes: Halfplanes in H axe filtered/replaced by their boundary lines by lo- 
calising the region containing the centre of min-disk either in the interior of 
these or in the interior of their complement. We do this by invoking LOCALISE 
with the set of boundary lines. In each iteration at leaist Vs of halfplanes 
are filtered/replaced in this way. 


Wedges: Wedges in W are filtered/replaced by localising the region containing the 
centre of min-disk in one of the four quadrants at the apex of a wedge (Fig. 20(d)) 
We then replace or filter the wedge, as the case might be. We call the four 
hnes, two sides and two normals, ^lssociated with each wedge to be mates of 
one another. 


I I. 


>'Uft 


iTrnrriB 
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We proceed exactly as in the case of line bisectors by first invoking LOCALISE 
with the set of derived lines of the wedges as input. Since for crossing derived 
line in the UR quadrant at most three derived lines of the corresponding 
wedge can intersect LL, we invoke LOCALISE twice. The first call is with the 
line/middle of these intersecting lines. Let LU be the region returned. We 
invoke LOCALISE again with the remaining line, if any, that intersects LU and 
whose mates do not. Let LL" be the region that is output. Thus we have a 
fraction of wedges in W whose associated lines do not intersect LLnLUfxLL" . 
This replaces/filters at least Vase of the wedges. 

We can get an 0{Nc log Nc) algorithm for computing the min-disk if we first convert 
all the wedges, halfplanes and rays to lines /points as above and then compute the 
min-disk of the resulting set consisting of lines and points only. To obtain a linear 
time algorithm, we do not differentiate between the bisectors of points and lines, 
associated lines of rays, halfplanes and wedges in the invocations of LOCALISE. 

We claissify the derived lines as active or inactive in a similar way as in the con- 
strained case. We assign weights to the active derived line in exeictly same way as 
the constrained case. At the beginning all the derived lines are active. We invoke 
LOCALISE with the set of active derived lines and if any active derived line does not 
intersect the localised region then we make it inactive for the rest of the computation. 
We replace/filter those objects in C none of whose derived lines are active. 

We repeat these steps until no update occurs. We later show that this adgorithm 
runs in linear time. 

We give a description of the algorithm for the unconstrained centre problem below. 
We will denote the active derived lines by V. 

Algorithm 3 UNCONSTRAINED-CENTRE 

Input: Bisectors of point -pairs and line-pairs, normals to rays, 
boundary lines of half planes, and normals to and sides of vedges 
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Output: Min-disk 
begin 

V < — input set 
do 

LOCALISE with V 
Update T> 

Filter/Replace objects in C 
while there is some update 
enddo 

Detenaine the centre by soae brute force aethod 
end I 


3.4 Analysis of the Algorithm 


In this section we establish that the algorithm runs in linrar tirw, First wr wij! 
analyse the constrained case and then the unconstrained one, 

In every iteration of the algorithm some of the wedges, halfplanrs. points and iina 
are dropped or some of the wedges, halfplanes and rays are replaced by points oj 
hnes. The DAG above (Fig. 22) illustrates this conversion. Further each localisatior 
of J with respect to some derived line reduces the weight by at least one fourth 
For example, when a wedge having four intersections with J is converted to anothei 
wedge with three intersections with J then the net reduction in weight is 16 
from 64 to 48. We show this reduction in Fig. 21. The number of intersections ar« 
denoted in the subscripts, and the number of objects resulting from the conversior 

IS wntten on top of arrows. Line pairs and point pairs are denoted by CV and VT 
respectively. 


Let Nu denote the number of objects of type U. 
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Figure 21: Reduction of Weights Per Iteration 
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Figure 22: Conversion eunong Object Types 
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Theorem 3.5 Algorithm CQKSTRAIMED-CEirTRE w Untar m Sr 4 ♦ V, ♦ + 

Nw. 

Proof: 

Let Wv be the total weight of active derived tinea. 

At least one of the lines derived from an object or a combinalitm ^^f object* in 
C intersects L. Otherwise we can straightaway filter or replace the i orrmjMinding 
objects. Then we have following bounds on Wp 

Wp > 27/Vw + 2?/V« + 21 Sn + 27.VJ2 4 l«.Vr . 

Wp < 48/Vvw + 21 Nn + 21 Nk + ISA'c + inSr. 

Therefore, 

W'n/48< Np + Nc + Sn -i- Nft + Nw < 2HV/27. 

Hence 

Wp « e(Np + Nc + /Vr + /V,^ + Nw)- 

In every iteration the derived active lines of at least half the total weight are made 
inactive. The new derived lines of at most three fourth of this weight are added. 
So at least a fraction (around one eighth) of total weight ta reduced. The pruning 
step is also linear in the number of objects which has the same complexity as the 
number of active derived lines. Therefore the algorithm is linear in Wp, Hence the 
algorithm CONSTRAINED-CEmE is linear in Np + Nc + Np + Nn ^ Nw m claimed. 
I 

The procedure LOCALISE, which locates the unconstrained centre with respect to 
constant fraction of weighted lines, uses the algorithm COHSTMIllEO-CEirmE as its 
basic routine. So, linearity of latter implies the linearity of former. We now prove a 
linear time bound for the algorithm inia)NSTRAIMED~C3nrRE. 



57 


Theorem 3.6 Algorithm UNCONSTRAINED-CENTRE is linear in Nw\- Nn Nji -\- 
Nc + Nv. 


Proof: 

Let Wd be the number of active derived lines in the algorithm UNCONSTRAINED-CENTRE. 
Similar to the constrained case above, at least one of the derived lines of each object 
or the combination of objects in C is active. Otherwise we filter or repl2ice the 
corresponding object/objects. Thus we have following bounds on W-o 

Wv > 27iVw + 27iV« + 27iV,? + 27iV£/2 + 18iV7>, 

Wv < 64iVw + 27iVw + 27iV^ + 18iV£ + 18 iV 7 ,/ 2 . 

Therefore, 

Wvl^< NT> + Nc^Nn^Nn-\-Nw < 2Wij/27, 

and hence, 

Wv = Q[NT> + Nc + Nn + NnA-Nv^). 

The pruning step is linear in the number of objects in C. Furthermore, in every 
iteration, at le<ist one eighth of the total weight of derived active lines is made 
inactive and new active lines of at most 3/4 weights are added, with the consequence 

that the total weight is reduced by 1/32- Therefore the algorithm is linear in VFp. 

Hence the algorithm CONSTRAINED-CENTRE is linear in iVw + N-h + N-ji + iV/: + N-p 
as claimed. | 


3.5 Concluding Remarks 


In this chapter, we have described an optimal algorithm for computing the smallest 
disk that intersects a finite collection of geometrical objects, containing convex 
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polygonal disks, wedges, baifplanes, liiw», poinU. Imr »rgtt»rnt» Ami rAV* In r<»m- 
bining the prune and search lechniqur of Mrgiddu w»lh ihr ninr-l iiJrA «.f frpjA* mg 
complex geometrical objects by simpler cmrs 

It would be worth investigating whether this approach ran Hr u*rti ui « otoputr the 
smallest intersecting disk of a collection of objects that mrlinirs simplr |Milygons 
and also look at higher dimensional generalisations of the proHlrtm studir^l hrrr 

In the next chapter we again see one more appliratiun uf the itsotHfirfl prune and 
search, which was used in this chapter, in the rornputatiun of a rrntrr f»r a 
finite planar set of points. 



Chapter 4 


Computing a Centrepoint of a 
Finite Planar Set of Points in 
Linear Time 


The notion of a centrepoint of a finite set of points in two and higher dimensions is 
a generalisation of the concept of the median of a set of reals. In this chapter, we 
present a linear time algorithm for computing a centrepoint of a set of n points in the 
plane, which is optimal compared to the 0(n log® n) complexity of the previously 
best known algorithm. We use suitable modifications of the ham-sandwich cut 
algorithm in [23] and the prune and search technique of Megiddo [21] to achieve this 
improvement. 


4.1 Introduction 


We all have an intuitive idea as to what phrases like “the very centre of the square” 
or “the very centre of the city” mean. To capture this intuition in a quantitative 
way, the centre of a set of n points, "P, in is defined as the maximal subset of 
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such that any closed halfspace intersecting this subset ronHuns at least \n/id + 1)] 
points of P [35]. This subset is non-empty for any finite ronfigiiratii.n of ikmiUs (see, 
for example, [13]). Furthermore, it is closed anti convex A crntrrpinni h a member 

of the centre of V. 

On the real line % a centrepoint is no other than a tnetlian of P. Thus a i entrepoint 
can be viewed as a generalisation of the median of a set of real.-* On the tiiher hand, 
the centre can also be viewed as a particular t-hull of P. The k-hull «f P is a maximal 

a 

subset (closed and convex) of such that any rloaed half spate inlersiecting this 
subset contains at least k points of V. For instance, the I hull of P i» its roovrx hull 
and the centre is its fn/(d+ l)l-huU. The properly of baJametl partititmiiig makes 
the centre point useful for efficient divide and conquer algorithm* in geximetrical 
computing and large scale scientific compuling[26, 32, 25, 3ti]. Hrt rntly Donoho 
and Gasko have suggested that centre point can be usetl aa “robust" and high 
“breakdown point” estimators for multivariate datasets [lOj. 

The interesting algorithmic problem of computing a rentrr{xiint has l>r<en considered 
by various researchers. Cole et al gave an 0(n log* n) algorithm for computing a 
centrepoint of a planar set of points [9]. Subaequenlly, {’«lr improvnl this bound 
to 0(nlog®n), using the powerful technique of slowing down a sorting network 
[8]. In this chapter, we propcae an optimal linear time algorithm for computing a 
centrepoint of a planar set of points by using suitable modiftcations of the ham- 
sandwich cut algorithm for a pair of separable point sets [23] and the prune and 
search technique of Megiddo [21]. 

Linear time algorithms, however, were known for computing an appwxtmote or c- 
centrepoint [19, 32, 23]. We obtain this weaker type of centrepoint if we decrease 
the lower bound, in the above definition of the centre, to fn{l - f)/(d + 1 )], where 
0 < £ < 1, Actually, Megiddo [23] only gave an algorithm for computing a partition 
of a (planar) set of n points with two lines such that each closed quadrant contains at 
least [n/4j points. An algorithm for computing an e-centrepoint, where 0 < £ < 1/4, 
is implicit in this. 
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The thesis proposes an optimal algorithm for computing a centrepoint of a planar 
set of points by using an interesting modification of Megiddo’s prune and search 
technique[21]. This consists of abiding a few extra points in each pruning step so that 
a subspace of the original solution space is retained, while ensuring a net deletion 
of points. In the description of our algorithm, we assume the usual RAM model of 
computation; the point set "P, however, is not assumed to be in general position. 

This chapter is organised as follows. In section 4.2, we discuss which points to prune. 
In section 4.3, we describe the method used to find these points. The algorithm is 
presented in section 4.4. Section 4.5 contains an analysis of the time complexity of 
the algorithm. Concluding remarks are given in section 4.6. 


4.2 What to Prune 


Let P be a finite set of points in the plane. In the subsequent discussion we use 
the following notations. We denote the centre of P by CEITIIE(P) and the A:-hull of 
P by HULL(fc, P). We use the notations Vh, 'PgH'> 'PfgHj •••to denote the points 
of P contained in if, G H if, F ClG f\ H,... respectively, where F, G, H, ... are 
any closed or open halfplanes. We denote the complement of a set S by As we 
frequently need to use the numbers flPl/S] and f|P|/3] — r|Pl/4] in this atnd the 
following sections, we denote these by N-p and Mp respectively. 

The basic idea of our algorithm is to use the prune and search strategy of Megiddo [21]. 
Clearly, we cannot hope to compute ceitre(P) by a naive application of this 
technique, since the centre of a reduced set need not be the same as the centre 
of the original set. However, it might be possible to prune points in such a way 
that the centre of the pruned set is a subset of the centre of the original set. If so, 
by repeated pruning we may at least be able to compute a centrepoint, if not some 
larger subset. Below we show that this is indeed possible, and as a first step towards 
this goal we make the following important observation. 
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Figure 23: CEITBEO? - T) C CEUlEtD 

Observation 4.1 IfT is the set of vertices of a triangle that contams CEfTtKl*). 
then ceitke(7^) is a subset o/ceitbb(T UP ). 

Proof: Let c be any centrepoint of V, i.e. c € CBitlE(P). By drfiiiitioii, any 

closed halfplane, say E, that contains c also contains at least Np points of P. Since 
by assumption c is contained in the triangle formed by T. we have T n // 0 so 

that H contains at least one iK>int of T. Thus H contains at least Np + 1 ~ A'put 
points of P U T. Since H is arbitrary, it follows that c is in clfTit{P U T). | 

The above observation has the following important consetiuenre. If we ran find a 
set of three points, T, in V such that the triangle formed by these points rontains 
the centre of P- T, then by discarding these three points we can achieve the goat of 
ensuring that the centre of the pruned set is a subset of the centre of the original net. 
The following lemma gives a sufficient characterisation of such a triplet of {>oints. 



Lemma 4.1 Let T 6e three points ofV such that miLL{Np -- 1,P) is contamed in 
the (closed) triangle formed by 7. Then ceitre(P - T) is a subset of CEITIUS(P). 

Proof: Let c be a centrepoint of P — T and T be the triangle foriiMHi by T. 

We claim that c lies inside T. Otherwise, if c lies outside 1\ and therefore outside 
hull(JVV — 1,P), then there exists an open halfplane that contains c and at the same 
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Figure 24: Removal of three points may expand the centre 


time contains less than N-p — 1 points of P. It can be easily seen that this halfplame 
contains less than Np — 1 = Np-r points oiV — T (Fig. 23). This contradicts the 
assumption that c is a centrepoint oiV — T. Hence c is contained in T eind therefore 
ceitre(‘P — T) is also contained in T. 

The proof of the result now follows from Observation 4.1. | 


Remark. We would like to point out a subtlety involved here. Had we chosen 
the triangle T to contain ceitre( 7^) instead, we could not have guaranteed the 
conclusion of the above lemma. Fig. 24 shows why. 

The above lemma suggests an algorithmic solution to the problem of computing a 
triplet of points that can be pruned. Since an open haJfplane that contains less 
than k points of V does not interesct BULL(A:,'P), we find three open halfplanes, eadi 
containing less than Np — 1 points of V and situated so that the intersection of their 
complements is a bounded trimgle. This triangle contains hull( 7 V 7 > — 1,V). If this 
triangle is of non-zero area then a required triplet is formed by choosing a point 
each from the closure of the pairwise interesctions of these halfplanes (Fig. 25). 

The snag in this solution is that there are configurations of points for which we 
cannot find such a triplet for any choice of these open halfplanes. An example of 
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Figure 25: Pruning of triplets T from V 



Figure 26: A pathological configuration 


sucb a configuration is shown in Fig. 26, where the points are evenly arranged on 
the circumference of a circle. 

To overcome this problem, we enlarge the scope of the above lemma, allowing for the 
choice of four points. For this we briefly review the concept of a Radm point Any 
set of at least four points m the plane can be partitioned into two disjoint subsets 
such that the intersection of their convex hulls is non-empty. A Radon point of this 
set is a point in this intersection. A Radon point of four points is unique when these 
points are vertices of a quadrilateral of non-zero area (Fig 27). 


Lemma 4.2 Let Q be any four points of 9 such that the (closed) convex hull of Q 
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Figure 27: Radon point(s) of four points p,^,r and s 


contains iim±{Np -1,V). Then ceitre((P - Q) U {^}) is a subset of ceitre(:P), 
where q is the Radon point of Q. 


Proof; Let c be a centrepoint of {V — Q)U {?}. Consider any closed halfplane 
H that contains c. Then, by definition, it contains at least ■^(7>-c)u{j} points of 
(P - (3) U {q). 

We claim that c lies in the convex hull, Q, of Q. Let p be a point that lies outside Q, 
and therefore outside HULl(JVV — 1,7^). It is then possible to find an open halfplane 
that contains p and contains less than N-p — l = iV( 7 >_e)u{,} points of (V—Q)V {?}. 
Hence p is not a centrepoint of {V — Q)\J {^}. Therefore c cannot lie outside Q. 

To complete the proof, we have to show that any closed halfplane H which contains 
c contains at least N-p points of V. Clearly, H contains at least iV(7>_Q)u{j} points 
of {P — Q)U {g}. Three different cases arise, depending on the relative positions of 
the points of Q. 


Case 1: The four points in Q form a non-convex quadrilateral. 

This case is a trivial application of Lenuna 4.1. The three convex vertices 
of Q form a triangle that encloses vmL,{Np — 1,P) and the concave 
vertex is q. Thus by Lemma 4.1 ceitre((P - fi) U {?}) C ceitre('P). 
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Figure 28: Substitution of Q by its lUdon point f 

Case 2: The four points of Q form a convex quadrilateral but their liadion point 
g does not belong to H. 

Since the quadrilateral Q and the faalfplane H both contain c, their 
intersection is non-empty. Thus H contains at least one of the vertices 
of Q and therefore at least Np = + 1 points of V. 

Case 3: The four points of Q form a convex quadrilateral and their Radon point 
q belongs to H (Fig. 28). 

In this case H contains at least two points of Q. We can therefore delete 
q from H and still claim that H contains at least Np points of P. | 


Thus in all cases H contains at least Np points of V. Since H is arbitrary, c is 
a centrepoint of ^ as well. Hence cbitbe(('P - Q) U {q}) is a subset of CEITM(P). 

I 


The above lemma is the cornerstone of our pruning mechanism. In the next section 
we will show how to use ham-sandwich cuts to malce a clever choice of four open 

halfplanes so that we can prune a fraction of the input set by repeatedly applying 
the last two lemmas. 
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4.3 How to prune 


In this section and afterwards, we use the words left, right, up and down, wherever 
these are unambiguous, to simplify the arguments. 

Suppose we choose four open halfplanes, call them L, U, R and D, (mnemonics for 
Left, Up, Right and Dovm respectively) such that each contains less that N-p - 1 
points of “P and its closure at least Np points, and together they enclose a non-zero 
bounded area. Why do we expect this choice to give us a triplet/quadruple of points 
Q satisfying the conditions of Lemma 4. 1/4. 2 ? We give an intuitive justification of 
this below. 

If the pairwise interactions of “adjacent” halfplanes ( i.e. L and 17, and R, etc.) 
are empty we would get a configuration as shown in Fig. 29. In this configuration, 
the total number of points in all the halfplanes taken together exceeds the total 
number of points in P by approximately one-third! This is impossible. So we 
might attempt to construct the four halfplanes in such a way that this excess 
is distributed evenly among the pairwise intersections of the adjacent halfplanes 
and thereby obtain approximately Mp triplets/quadruples of points satisfying the 
conditions of Lemma 4. 1/4.2. 

It is possible to do this as the construction below shows. 


4.3.1 Computation of Open Halfplane L 

We fix Ir as follows. We determine an extreme point poiP with minimum abscissa, 
and join all the remaining points to it. We compute the line that passes through 
p such that its slope is the Np — 1-th largest of the slopes of the above lines. The 
open halfplane above this line is chosen to be L. Clearly, it takes linear time to 
compute L. This way we make sure that L contains less than Np — 1 points and its 
closure contains at least these many points of P. Moreover, this haJfplane contains 
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Figure 29: The intuition behind pruning 


at least two points on its boundary. Really speaking, this latter requirement is not 
necessary but it helps us to treat the four haJfplanes uniformly in the analysis of the 
algorithm. 


4.3.2 Computation of Open Halfplane U 

Since the point set "P can be degenerate we need to be careful in the construction 
of U in order that none of the closed quadrants determined by the boundlaries of L 
and U contains too few points of P. This will also result in an even distribution of 
points among the pairwise intersections of adjacent halfplanes. 

To achieve this we use the ham-sandwich cut algorithm of Megiddo [23]. However, 
the ratios in which we propose to divide the point sets are arbitrajry. As we show 
below, Megiddo ’s algorithm can be easily adapted to take care of this aspect. 

As is usual, we consider the dual problem, letting the bondary of I* be the y-axis in 
the primal plane. Under the duality transformation that we consider, points that 
lie on the boundary of L map to horizontal lines; the set of points that lie in L 
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U 


Figure 30: Computation of L eund U 


map to a set of lines, with negative slopes; ajad the set of points that lie in 
interior of L map to a set of lines, with positive slopes. In order to get the type 
of ham>sandwich cut we want, we put the horizontal lines in Let us assume 
that jC+ contains at least Np — Mp lines in the dual plane. It can be seen that the 
boundary of U corr^ponds to that point in the dual plane which hais less than Mp 
lines of £“ above it; at least these many lines of C~ passing through or above it; 
less than Np — Mp lines of above it; and at least these many lines of passing 
through or above it. We compute this point as follows. 

It can be easily seen that we can resolve a query for any positive (negative) slope 
query line in Megiddo’s method as follows. We first compute the Afp-th {Np — Mp- 
th) intersection of the lines in C~ (£"*■) with the query line. We then count the 
number of lines in {C~) lying strictly above and the number of lines in (£“) 
passing through this intersection point. If the sum of these two counts is smaller 
than Np — Mp (Mp) then the solution point lies below the query line. If the first 
count is greater than or equal to Np—Mp (Mp) then it is above the query line, else it 
is same as the intersection point. By thus dianging the method of query resolution, 
we get the solution point in linear time. This gives us the required halfplane U in 



70 


the primal plane (Fig. 30). 


4.3.3 Computation of Open Halfplanes D and R 


The open halfplane D is determined with respect to L in a similar manner by 
changing “above” to “below” throughout the above discussion. Thus we make sure 
that D contains the “lower end” of the boundary of L whereas (• contains the “upper 
end”. 

The halfplane R is also determined similarly except that U plays the role of L here. 
We also ensure that R contains the “right end” of the boundary of £/ whereas L 
contains the “left end”. 

The idea behind choosing D and R in this way is to make the boundaries of the 
halfplanes U, R, D and L the adjacent sides of a bounded quadrilateral such that 
the halfplanes face “outwards”. This is so because if the interior of I contains at 
least N-p — Mp points of "P then the boundaries of both U and D are disjoint from 
the boundary of L and it can be seen that the intersection of the complements of 
the above halfplanes is bounded. As explained before, it can also be seen that jRn D 
contains non-zero points of V. Moreover, since each of these half planes contadn 

less than Np- I points of V, the intersection of the complements of these encloses 
wsii.{Np-l,V). 


However, in the computation of U it is quite possible that the set £+ contains less 
than Np Mp lines. The consequence of this is that the computed U may have 
the same boundary as that of L and thus the intersection of the complements of the 
open half planes may be unbounded. This is inadmissible in our algorithm. So we 
need to take care of this degenerate case separately. 
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4.3.4 The Degenerate Case 

If there axe more than N-p points on the boundary of L then the boundaries of U 
and Z), as computed above, is the same ais that of L. This degeneracy is detected in 
the algorithm when £+ contains less than N-p — Mp lines i.e. interior of L contains 
l^s that Np — Mp points of V. There is no loss of generality if we assume that L 
also contains less than Np — Mp. Otherwise, we can switch the sides of L and let the 
interior of X be our new L. We will then be able to compute the other halfplanes, 
as required, with respect to this open half plane. 

Let the interior of X be the open halfplane R. We compute the open halfplanes U 
and D as follows. We first distribute the alternate points of 'P on the boundary of 
L between the sets R and L. Let these resulting sets be Sr and Sl respectively. We 
then compute U and D by the ham-sandwich cut algorithm such that these satisfy 
the following properties: these contain less than total Np — 1 points of P and contain 
am equal number of points of Sr and Si,. The halfplanes U and D are computed such 
that these contain the “opposite ends” of the boundary of L. It can be seen that if 
X, jR, U and D are computed in this manner then each of the pairwise intersections 
of the adjacent halfplanes contauns at least Mp (approximately 2Mp) points of P. 
Since the intersection of the complements of these contains HULL(lVp — AZp,P), we 
cam apply pruning at least Mp times similarly as in the non-degenerate case that 
we describe below. 

The degeneracy of open halfplanes taken caxe of, we may assume safely that the open 
halfplanes X, £/, R and D can be computed such that these meet our requirements. 


4.3.5 The Pruning Step 


We have been able to ensure by the construction of the halfplanes as above that the 
closure of each of the sets LOU, Ln D a.nd ROU contains at least Mp points of 
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V. We shall prove later that the closure of /2n also contains at least these many 
points. However, for the rest of this section we will assume this. 

It is now clear how we can prune points. Two points of detail, howevrr, must be 
noted. First, in order to ensure that the conditions of Lemmas 4.1 and 4.2 remain 
valid throughout the pruning step, we must choose a triple or a quadruple of points 
in such a way that, whenever there is a possibility that the conditions of the above 
Lemmas are violated in the successive pruning steps, we delete an interior i>oint in 
an open halfplane. Second, to maximise the number of points that are pruned, we 
must ensure that no two points of a quadruple, selected for pruning, belong to either 
Rf\LoxUr\D. 

To implement the above observations we maintain the points that are candidates for 
pruning in six disjoint sets, viz., Vjju, VuRt 'Rrd-, "Plr «*^d Pyo. The points 
on the boundaries are put in the relevant sets. So, the four sets, Lf\l/ , LClD, ROU 
and RnDf are now effectively divided into six sets, three of which correspond to 
choices of triangles T. 

We discard the triangles T and substitute quadruples Q by their Radon points 
sudh that a maximum number of the above half planes contain an interior point. 
Substitution of Q is done as follows. If Q form a convex quadrilateral we delete it 
from P and add the intersection point of the diagonals to P. Otherwise we delete the 
convex vertices but retain the concave one. We can repeat this pruning procedure on 
the reduced set of points thus obtained, since the halfplanes L, £/, D and R continue 
to contain less than N-p — 1 points of the reduced set P, until one of the four sets is 
empty. We note that this reduces the size of P by approximately one fourth. 


4.4 The Centrepoint Algorithm 


It is now clear from the discussion in the previous sections, how we can find a 
centrepoint of P. 
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In each iteration we compute the points that are to be discarded or replaced. By 
throwing away these points we reduce the size of the set by a non-zero fraction. 
When the size of the set becomes so small (of size at most 10) that no more points 
can be discarded we halt the pruning procedure and compute a centrepoint by smy 
straight forward method. 

The algorithm is given below. 

Algorithm 4 CENTREPOINT 

Input: Set of Points V 
begin 
do 

Compute the open halfplanes L, U, D and R 
Update V by deleting T and replacing Q by 
their respective Radon points 
while there is some replacement/deletion 
enddo 

Compute a centrepoint by any bmteforce method 
end I 


We justify that anytime during the pruning step each halfplane contains less than 
N-jt-l points where R is the current set of points. 

We first argue for a deleted point that is also an interior point of a haJfplane. We 
first consider the case when four points axe pruned and their Radon point is added 
to the set T. If a point that is pruned lies in the interior of a halfplane and the 
Radon point does not lie in that halfplane then the number of points in this half 
plane is decreased by one. Since the total number of points decreases by three, the 
number of points is less than - 1 = - 2 where H is the new set. Now suppose 

that the Radon point also lies in this halfplane. Since a halfplane that contains a 
Radon point of four points contains at least two of these points, in this case also the 
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Figure 31: The woret ca*e of pruning' 


number of points in the halfplane decreases by one. The same argument holds for 
triplets of points. 

Now we consider the case when there is no interior point among the four pruned 
points in a halfplane. In this case the Radon point aJso does not belong to it. We 
consider the worst case in which for every choice of a quadruple or a triplet of points 
those chosen from this halfplane lie on its boundary. This case is simpler to analyse 
and it is easy to see that same argument extends to the other cases. H is clear that 
when we start pruning, the number of points in each of the open halfplanes is less 
than N-p — 1 points. Let us see what happens when we have to prune the last triplet 
or quadruple of points after applying pruning Mp ~ 1 times. The maximum number 
of points in this open halfplane is less than 


Np — Mp 


since Mp points in the corner quadrants are already pruned. 


The total number of points at this instant is 

\V\ ^ZMp + Z 


The number in the previous expression is exactly one less than the ceiling of this 
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number. Thus the last set of points can also be pruned. It can easily be seen that 
in the intermediate steps also we can do the pruning. 

The proof of correctness of this algorithm now follows from Lemmas 4.1 and 4.2. In 
the next section we give an analysis of the running time of this algorithm. 


4.5 Analysis of the Centrepoint Algorithm 


For the purpose of the proofs in this section, we assume that the boundaries of the 
half planes L^U^D and R do not contain any point of V. This assumption is not 
necessary to establish the linearity of our algorithm but it simplifies the proof to a 
great extent. This can be achieved by slightly perturbing the points of V such that: 
no point migrates across the boundary of any halfplane; each corner region, such as 
Lf\U etc., contains exactly M-p points of P; and each halfplane contains at least 
Np points of V. Intuitively, such a perturbation does not matter because if we 
can prove that the perturbed set ViiD contains Mp points then these many points 
also belong to the closure of P/uj in the non-perturbed set 'P . As a consequence 
we do not need to state explicitly whether the halfplanes are open or closed in the 
following discussion. 


To prove that the algorithm is linear we have to show that the size of V is reduced by 
at least a fraction in each pruning step. We know by construction that each of the 
sets, Vwi 'Pur Juad Pld-, contains Mp points. We have to show that Vrd also con- 
tains at least Mp points in order to prune at least these many triangles/quadruples 
from V. For this we will have to consider all the possible relative positions of the 
halfplanes L, f/, R and D. Several cases arise. A few of these can be straight away 
discarded by using the following lemmas. 


Lemma 4.3 Let F, G and H be three halfplanes. Then F CiG r\ H ts a bounded 
triangle if and only if F OG C. H- 
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Figure 32: Interaection of three halfphuiee 


Proof: Straight forward(Fig. 32). I 


Corollary 4.1 Fn^CiH is a hounded triangle if and only if F C\G C W if ond 
only if F f\ H cG. 

The sets 17, L, R and D satisfy some additional constriunts also on account of their 
specific methods of construction. 

Let us denote by pan the intersection point of the boundaries of any two halfplanes 
G and H. For the sake of simplicity, we may assume without any loss of generality 
that no three boundaries of the above halfplanes intersect at a point. If they do, 
then we can treat the said configuration in one of the cases discussed later on. 


Lemma 4.4 The intersection of the halfplanes U and D is contained in L if and 
only if Pud is contained in L, i.e. U f\ D C L 4=» PUD^ 1“ 

Similarly, LORcU pm € U. 


Proof: We prove only the first equivalence. The proof of the second is similar. 
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(=») Easy. 

(<=) Ixjt Pud be in L. Then four cases arise, depending on the orientation of the 
half planes U and D. 

Case 1: V dV C L 
Case 2; UnDcL 

These cases are not possible. For by construction, the halfplanes U and 
D contain the “opposite ends” of the boundary of L. 

Case 3: V nV C L 

This case is also impossible since it implies that the centre is empty. 
Case 4: U f\D C L 

This is the only permissible case{Fig. 33). I 


Hence proved. | 

We have the following similar lemma when puD H®® 

Lemma 4.5 The intersection of the kalfplcnes U and D is contained in L if and 
only if Pud contained in L, i.e. U f\D C L < — > puD € L. 

Similarly, LC\ RcU Plr € U. 

Proof: Here too we prove only the first equivalence. 

{=») Easy. 

Let PUD be in L. Then four different cases arise. 
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Figure 33; Either U r\DC L ot U n D cl 


Case 1 : U ClD C L 
Case 2 ; llnD cT 

These cases are not allowed by construction. 

Case 3 : V C\~D clj 

To see that this case is aiso not permissible, let us count the number of 
points in Vi. First we prove that Vi =s u Vm. 

= 'Plud U ^£^25 U Viffoy 
= Vw U ViD^ 

since Vlud U = Vld and VwD U 'P^yu = Vw Hence, 

\Pl\ = iPu/UPiol 

= I'PuLfl + \'Plo\ - {VijudI 
— M-p + Mp IVijud] 

< 2Mp. 
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This contradicts the fact that Vi contains at least N-p points. Thus this 
case is also not possible. 


Case 4: V Ci D C.X 

This is the only permissible case. | 


Hence proved(Fig. 33). | 

The above lemmas have the consequence that either UC\DcLoxUC\D(lL, and 
similarly either LCiRcU at LOR CU. Now we can prove the following theorem. 


Theorem 4.6 There are at least M-p points in Vrd- 

Proof: For the proof, we again consider all the possible relative positions of the 

haifplanes £, U, R and D. The following cases arise, depending on which of the 
four quadrants formed by the boundaries of U and L contains Prd' 


Case 1: pm ^ UHL (Fig. 34) 

In this case £> n i/ n i and i? n n L are non empty. Therefore, from 
Lemma 4.4, 

UnDcL and LORcU. 

Using Corollary 4.1 these respectively imply that 
U ORcT and LOD CU. 

Let X be a point oiROD. Then x lies in exactly one of the sets, UOL, 
U n L or Ur\L. Since VoRcL therefore x^UOL. Similarly 

X ^Unl. Now ROD is convex and there is a point pm that lies in 
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U C\ L, therefore x does not lie in 17 n Z either because a convex region 
cannot intersect only the opposite quadrants of a pair of lines. Thus 

RnDcUDL, 

and hence, by Corollary 4.1, 

17 n jR c 75 wid X n £? c 1?. 

Now, 

V = 

= U U U Ppjyjj 

= 'Pl U 'Piu U U 'Pijjji since UCiRcX 

= “Pl U Viu U U Pjr^, U PixulTS since X H £> C 77 and Tr\D 
As these axe disjoint sets 

\P\ > {Np) + {N'p — M'p) + (Np — Mtp) + {Np — Mp) ’i- I'Pijyjijfl 
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> rm/3i+3rwi+i^i7S2?i 

> 1 ^ 1 , 

which is a contrewiiction. Hence this case does not occur. 

Case 2: prd qT7 C\ L (Fig. 35) 

In this case/inLnI7^0 therefore ROLcU (Lemma 4.5). We will 
have to consider four subcases. 

1. RnDcL 

2. i 

Since RC\D <Z L therefore R and D contain the “opposite ends of 
the boundary of L Since RdLcV therefore from Corollary 4.1 
H and 1/ contain the “opposite ends” of the boundary of L. This 
implies that D and U contain the “same ends” of the boundary of 
L, which is not possible by the construction of D. 

3. 'MnDcL 



By Lcnnni& 4.3 this implies. 


InDcR. 


Therefore, 


Thus 


pRD - 'pRDL^'pRvil 

\'Prd\ '> Pol! 

> /Vj> — Mf>- 


4. RC\'BcL 
This implies 


RC\VCD. 


Otherwise consider a point * € /? O (/ O 25, 


x^Rr\Ur^lD ^ ^ x£ L 


Again, 

x€jR =► X ^ Rf\L s» X E.1J. 

which is a contradiction. 

Therefore, 


Hence, 


'Prd = Prdu U Prqjj 
- PruUPrdu- 


\Prd\ > I'pRul 

> Mp. 


Since prd E L so these four cases exhaust all the possibilities. 
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Figure 36: prd G U r)L 


Case 3: piu> € U nX 

This case is similar to the previous one. 

Case 4: pno € fl X (Fig. 36) 

Again four different subcases are possible. 

L Rn L CU mdUnD C L 

Consider a point x € Rf) D. If x € L then x € U hy ROL cU. fi 
X € I then xeUhyUnDcL. Since ROD is convex, it cannot 
intersect exactly two opposite quadrants of a pair of lines. As there 
is a point pmi in 1/ n L hence Rf) D (Z U C) L, 

Assume that IViwl = Then, 

'Put = 

= :%«dI-'^SdU7Vj®'J%ibb, since £7 nilc land 
= T’jicUPtED'J'PtijaiUPprSB. since JJn Pet/ 

As the sets on the left of the last equality axe disjoint sets, 

IPutl ~ + I'PekdI + 



Substituting the values of the different terms, we get, 

|'P| — 2JVp + M 7 > > m + {N-p - Mp - m) + {Np - Mp - m) 

I'PulrdI ^ m+\V\-4Np + ZMp 

< m - Mp - + \Vl 

Since the size of the set is non-negative, 

m > Mp. 

2. ROLcU mdUnD cl 

Since pnn eUnZ therefore pm, lies in U. There are four different 
possibilities. 

• RCxDcU 

• "RCxV C U 

By a similar argument as in case 2 we can show that the above 
two cases are not possible. 

• RnDcXf 

Consider a point x C Rf)L. 

RDLCU => xCU 
Rn'D cW xeD 
UnDcI => x€l 

X € Ln L, 

which is a contradiction. Hence RdL = 0. Therefore Rr\L C TJ 
and U Ci D C L — a subcase we shall consider later. 

• RO DcTJ 

Let IVrz}] = m. We compute 
Now, 

V = 
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= 'Pd ^ 'PpR U "PoSi; ^ P'Dl^w 

= Pd U P 5 ;j U Pgj; U PoKiTi U V^Ruz, since RDU CD 
~ "Pdr ^ 'Pm 'Pdul ^ P'd'rUli 


since Lr\U C R. 

Substituting the values of different terms in the above equation 
we get 

|P| ^ — m + N-p — M-p + Np — 2Mp + I'P^rulI 

\PBJiTJz\ ^ ~ + |P1 

< m-M^-4nPl/4l + IPl. 


Since the size of the set 'P^rul negative therefore, 

m > Mp. 

3. Rf\LcV !iXi^Uf\DcL 

This case can be dealt with in a similar way as above. 

4. ROLcV^d UDDcI 

Since pro € U here also we shall consider four different cases. 

• RdDcTI 

• RoVcTJ 

These cases do not occur for the same reasons as discussed 
above. 

• RnDcU 

Prom Corollary 4.1, 

Rnu c i?n J?. 

Since RDU contains Mp points of P, 

|PrdI > \pRu\ 

> Mp 



• ^nDcU 

We need only consider the case Rn'D cT since other cases, 
viz. ROD cT,'Rf\DcT and 7? fl 15 C I, can be dealt 
symmetrically as above (actually they leads to a contradiction 
when considered along with 71 0 Z? C U)- Let Prd contain m 
points. As before we compute the number of points in Pjj^ux- 

P = Pd^Pu 

= PdUPsrUP^^ 

= Po'JPsR^Ps-mj^'hsJiU 

= Pd Pj5r U Pxiu 'PshUl *-> '^TiWl> ^ fl £/ C 75 

since LC\D <zR. 

Substituting the values of different terms in the above equation 
we get 

\P\ ^ ZVp + N-p — m H- Np — Mp + t^p — 2Mp + 

^ + 3iWp + I’Pj 

< m-Mp-i\\P\l4:\^\P\. 

Since the size of the set Px^xHH non negative therefore, 

m > Mp. 

Thus theorem is proved for the last case prd € 77 H X also. | 


nee proved. | 


mbining the earlier theorems and lemmas we get the following result. 
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Theorem 4.7 A point in the centre of a set can be computed in linear time. 

Proof: In each iteration at least ZMp (~ 17^1/4) points are deleted. 

If T{n) is the running time of the algorithm for an input set of size n (IPI = n), 
then it satisfies the following recurrence. 

r(n) < max T{n — 3k) + 0{n) 

kyM'p 

=» T(n) < f(n - 3Afp) + 0(n) 

Since T{n) = 0{n) from the above recurrence relation, the claim of the theorem 
follows. I 


4.6 Concluding Remarks 

We have presented an optimal algorithm for computing a centrepoint of a finite set 
of points in the plane, thus providing one more example of the power and versatility 
of the prune and search paradigm. 

It would be worth exploring how tUs speeds up algorithms which uses the centre- 
point computation as a baisic subroutine. 



Chapter 5 


Designing Algorithms Using 
Partial Sorting Networks 


In this diapter we present a general technique that is a sequel to the parametric 
searching technique of Megiddo. The latter technique is used to design serial 
algorithms with the help of efficient parallel ones. However, these algorithms are 
not very efficient for the class of problems that use parallel sorting networks in their 
solutions. We modify this technique so that we obtain optimal linear time algorithms 
for some of these problems. We do this by synthesising prune and search technique 
and parametric search technique. 


5.1 Introduction 


The parametric searching of Megiddo is as follows — Let there be an efficient parallel 
algorithm for a problem A such that solution of A can be used in the solution of 
another problem S. Then, in some cases, we get an efficient sequential algorithm 
for JB by exploiting the efficient parallel mechanism of the parallel algorithm for 
A. This technique has been applied to a wide variety of problems yielding efficient 
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s-igorithms. In particular, we achieve good results for the parameterised problems 
that use parallel sorting algorithms in their solutions. We do this by replacing 
the evaluation of the parallel comparisons of an iteration in the parallel version by 
simultaneous resolution of these in the serial version. The running time of these 
algorithms is further improved by introduction of weights in the comparisons [8]. 
We then simultaneously evaluate at least a fraction of total weight of these in every 
iteration to design more efficient algorithms. 

We further improve the running time of the above algorithms by introducing, wher- 
ever applicable, prune and search in these. We do this by seeking to compute the A:-th 
largest element of the input set instead of seeking to sort it as in previous techniques. 
For this, we run the sorting algorithm on a given input for a few iterations and then 
prune the set. It can be eeisily seen that this approadi is useful only where pruning 
is applicable. 

This chapter is organised as follows. In section 5.2 we review a few definitions and 
concepts related to the AKS sorting networks. In section 5.3 we present a method 
of applying prune and search using these networks. In section 5.4 we optimally solve 
the problem that was first posed by Megiddo. In subsequent sections we apply this 
technique to linear programming problem, computation of ham sandwich cuts and 
computation of centre points in d-dimensions. 


5.2 Preliminaries and Definitions 


The AKS sorting network sorts a sequence of reals in 0{n log n) comparisons such 
that in each parallel iteration 0{n) comparisons can be performed so that there 
are in total 0(log n) iterations. A comparison operation in the network consists of 
comparing the contents of two registers and interchanging these if their specified 
order does not match. 

We can view the process of sorting in the AKS sorting network as the movement of 



gisters in a complete binary tree of n leaves. The order of the leaf (topmost) nodes 
this tree corresponds to the ordered sequence of n numbers. All the registers are 
itially in the root (bottommost) node of the tree and in every iteration, most of 
e registers in a node move up whereas only a small number of registers move down, 
hen each of the registers reaches a leaf such that each leaf has only one register 
esent in it then the computation is smd to be over and the elements in the leaves 
e in sorted order. 

e use the same notation as in [1] by Ajtai et al. Let be a set of n registers, 
aich contain elements of the input set, in the AKS sorting network. Let the set 
registers assigned to a node t of the complete binary tree in a-th iteration be 
inoted by S^{t). Thus, S°(<) is if t is the root of tree, otherwise it is 0 (empty 
t). The restriction of the function 5® to a fixed level of height h is denoted by 
*’*. 5® is defined in such a way that each of the nodes at any height, A, contain 
me number of registers. Let N{S°'^) be this number. We bound this number in 
e following theorem. We reassign the constants qi, qj and Cj used in the paper {1] 
q, Q and c respectively in the following discussion (Q < f < I/c < 1). 


heorem 6.1 ([l]) The number of registers in cny node at height h in a~th itera- 
on has the following properties: 

1 . 

2. i/?®-'^2-®-*n < c then N{S^'^) = 0 otherwise > 9 «~* 2 -®"*n. 

'e can see from the above theorem that the number of registers in any node in 
•th iteration is bounded, both above and below, by terms that are in decreasing 
sometric sequence. The AKS sorting network has another useful property. The 
dering of the elements after each iteration is an approximation of the sorted order 
the input set. Moreover, at any instant the number of the elements having a 
rge displacement is very small. We bound this number in the following theorem. 
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Let the relative position of a number a; in a sequence be defined as the fraction of 
numbers in the sequence that are sm<dler than x. We also assume that the positions 
to right correspond to the laxger elements. 


Theorem 5.2 ([1]) Let G^{t), where i is a node at height h (h < a), be the set 
of roisters in the nodes, right of t, at height h such that the relative positions of 
elements in these are smaller by at least p than the relative position oft. Similarly, 
let Hf,{t) be the corresponding set of registers with larger elements on the left of t 
such that the relative positions of the elements contained in them are larger by p 
than the relative position of t. Then, 

and, similarly, 

where M = and m > 1. 


Now, we seek to apply the Megiddo’s technique to optimisation problems with the aid 
of these theorems. In these, we frequently need to prune redundant hyperplanes in 
E^. The multi dimensional prune and search procedure proposed by Megiddo is used 
for this purpose [22]. According to this, we can prune a fraction of the hyperplanes 
by locating the optimal point with respect to a finite set of query hyperplanes. The 
following theorem establishes the existence of these query hyperplanes. 

Theorem 5.3 ([22]) For any dimension d there const constants A = A{d) and 
B = B{d), with 0 < B< sudi that A queries suffices to determine the position 
of optimal point relative to at least Bn of n hyperplanes in E'^, where 

A{d) = 2^-^ 


and 


B(i) = 2*-"'. 
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It is clear that this procedure of pruning hyperplanes is applicable only if there 
exists an oracle, say that takes linear time to determine the relative position of 
the optimal point with respect to a query hyperplane. 

In the next section we discuss the pruning of elements in the input set. 


5.3 Pruning Using the AKS Sorting Networks 


Let 5' is a set of n elements. Then, there aje altogether N (= ’Cj) distinct 
comparisons possible between pairs of the elements of S. We assume that the 
output of any comparison depends on the value of a vector x i.e, it depends on 
the membership of x among disjoint subsets of the domain. In the problems we 
discuss, these are the regions determined by hyperplanes. We call these critical 
hyperplanes. If the optimal points x" is known then the complete order of elements 
of S is known. The critical hyperplane of each unresolved comparison intersects the 
localised region. Furthermore, it is to be noted that localisation of x“ with respect 
to one of these hyperplanes as well as localisation with respect to 0(n) of these 
(only a fraction of these will be resolved at a time using multi dimensional prune 
and search) require linear time( Theorem 5.3). 

We pose a hypothetical problem of computing the k-th largest element, s*, of S at 
X . What are the elements that can be pruned ? Clearly, these are the elements 
that are not the ^-th largest element at x*. We formulate our problem as follows: 


Subproblem 1 How can we determine, in linear time, 0{n) elements o/$ that are 
not same as s* given an oracle ft that takes 0(n) time to answer the following query 
~ which side of a given hyperplane does x* lies? 


The real crux of the problem now is to identify the elements that have large dis- 
placement from Sk. For determination of these, we generalise the method of Lo and 
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Steiger in [18]. We determine a prismoid, whose base is the localised region, about 
the fc-th largest element such that most of the elements do not intersect it. The 
elements that do not intersect this prismoid are always at a non zero finite relative 
displacement from Sjt in the localised region and therefore, can be pruned. 


5.3.1 What to Prune ? 

Let <S be a set of hyperplanes in E'^. Suppose we can localise x* in a region J such 
that the total number of unresolved comparisons is only a small fraction of JV. We 
construct an infinite prism with J as its base such that only few of the intersections 
of any pair of elements of S intersect J. Then we prune the hyperplanes by the 
following theorem. 

Theorem 5.4 Let V be the prismoid determined by k — {d — !)[£«] tind k-\- {d — 
l)fcn] laryest intercept on each vertex of localised region (which is a simplex). If 
there are less than unresolved comparisons in the localised region then toe 

can prune approximately (u — elements in the computation ofk-th largest 
element. 

Proof: unresolved comparison corresponds to the intersection of a pair of 

hyperplanes in 5 which also intersects the infinite prism of which P is a part. 

We compute the minimum munber of these intersections needed so that the surface 
formed by ifc-th largest elements (fc-level) may lie above or below the prismoid at an 
interior point p. We only analyse the case when the surface lies above V at p. ' 

Let us consider the following worst case configuration. In this configuration, eadi of 
the hyperplanes that is above the fc-level and below the A: + (d - 1) [en] -level at any 
of the vertices Vf, with 0 < t < d of the localised region is above the prismoid at p. 
Let the set of these hyperplanes be H. Also, in the minimal case, no hyperplane that 
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is below the Jk-level at any of the e/s is above the prismoid at any point. Otherwise, 
the number of intersections in J increases. Similarly, no hyperplane, that is not in 
H but which is above the it + (d - l)f£nl -level at any of the v,\ intersects prismoid. 
Furthermore, none of these hyperplanes intersects any other hyperplane over J. 

Now, the hyperplanes in H in this configuration can be divided into several equiv- 
alence classes defined as follows. If two hyperplanes intersect or do not intersect 
the prismoid at the same vertices then these belong to same equivalence class. 
Each equivalence can be represented by a set of vertices of the base, at which the 
hyperplanes of this class intersect V. We call these as inttr$tcting vertiets of an 
equivalence class. A hyperplane in one equivalence class intersects a hyperplane of 
another equivalence class if their intersecting vertices are not related by the relations 
of either subset or superset. , 

If there are d equivalence classes of size fen] such that each has only one intersecting 
vertex then the number of intersections above J is minimum. The count the number 
of intersections in this configuration is If the number of intersections is 

less than this then the Jfc-level will lie inside the prismoid. 

We can prune those hyperplanes which are above or below the prismoid at all the 
vertices v,-. The maximum number of hyperplanes that will intersect the prismoid is 
when the equivalence classes are chosen as above and the intersections are divided 
equally between the hyperplanes both above and below the ib-levcl. Hence, we can 
prune the rest of approximately (n - d[-s/2en \ ) hyperplanes in S. | 


We can prove another similzu: theorem. 


Theorem 5.5 If there are less than ‘fc72fen]^ unresolved comparisons in the feasible 
region then the prismoid is formed by the largest (smallest) and {d- largest 

(smallest) intercept on each vertex of feasible r^ion (which is a simplex) and we can 
prune (n — dfen]) elements in the computation of the largest (smallest) element. 


95 


To use prismoid method we need to ensure that number of unresolved compaxisons in 
the feasible region should only be a small fraction of N. Then we can guaxuntee the 
existence of a non zero fraction of prunable elements. In next section, we compute 
the minimum number of iterations of AKS sorting network needed to ensure this. 


5.3.2 Why can we Prune ? 


We use the properties of AKS sorting network stated in the previous section to 
prove that the number of unresolved compaxisons in o-th iteration is bounded by a 
number that decreases in a geometrical progression with a. 


Theorem 5.6 The number of unresolved comparisons in the AKS sorting network 
in the a-th iteration, out of a maximum of "C^ ones, is bounded from above by 




Proof: Let us first count the number of unresolved comparisons of an element x 

in a register present in a node at height h in o-th iteration. Let p, be 2/c and t be 
that node which has the same relative position as x. Then 

2"‘/M < li, 

2 ” < iiM 

< 2M/c 

< 29 “"'^ 2 -“"^n/c 

< g“-*2"“n/c. 


Since > c, therefore, m = 1 satisfies the above inequality. 

Thus, the number of unresolved compaxisons of elanents at height k with a larger 
relative displacement than p from x is bounded by 
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This quantity when summed over each level of the computation tree is at most 
, 0<k<a 

< 2(-i-)2-"nl3. 

1 - q 


The number of remaining unresolved comparisons of x with elements having relative 
displacement smaller than ft is at most 2/u = Anfc. We take half of this in the total 
as each of these unresolved comparisons is counted twice in the sum. Further, we 
also add unresolved comparisons corresponding to any two registers in the same 
node to the computed total. Therefore, number of total unresolved comparisons is 
at most 


f(-i-)2-“n<2+2n/c'ln+ V 2^N(S‘’-^f/2 

/ Q<h<Oc 

< + 2n/c) n + 


Hence proved. | 


We can make this fraction as small as we need by choosing ot, q, Q, c appropriately. 
We can thus bound the number of unresolved comparisons in each iteration. But 
there still remains an algorithmic problem. Resolving all comparisons of even 
only a first few iterations of the AKS sorting network requires O(log n) steps of 
Megiddo’s multi dimensional prune and search. This will take at least O(nlogn) 
time. However, this is not agreeable. So, in order to make the pruning procedure 
linear, we have to leave another small fraction, say 5/2, of comparisons unresolved. 
We observe that this adds at most 2{Sn/2)n in the above quantity i.e. two elements 
for each unresolved comparison (and these may remain uncompared with all the 
other elements). Then, the total number of unresolved comparisons is 



^) 2 -» + | + 5 )„» 
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where S is also chosen appropriately. Now, we can apply Megiddo’s multi dimen- 
sional prune and search on the critical hyperplanes until only 6n/2 comparisons of 
the total comparisons of first a iterations are left unresolved. Clearly, this takes 
linear time. 

The method of application of the technique will become more clezir in section 5.4 
when we apply our technique to a few problems. 


5.3.3 How to Prune ? 

We can view AKS sorting network as having O(logn) iterations with at most n/2 
comparators in every stage. We give a weight of 4“'^ to each comparator at a depth 
j as in [8]. We choose the values of q, Q, c, a and £ appropriately and assign 
weights to comparators of first ct stages only. The construction of first a stages of 
AKS sorting network takes linear time. We give details of the algorithm to compute 
x“ below. 

Algorithm 5 Prune And Search, on AKS Sorting Network 

Input: S 
begin 
do 

Let C be the set of unresolved comparisons 
and W its total weight 
do 

Solve (d — l)-A(d-l) queries using ( Theorem 5.3) 
Resolve at least B{d-l) comparisons 
Update C and W appropriately 
while |(7| > 5n/2 
enddo 

Prune the hyperplanes in 5 ( Theorem 5.4) 
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Update k aad S accordingly 
while S is modified 
enddo 

Apply any brute force method to compute x 
end I 

In the above algorithm (d - l)A{d - 1) queries are put to the oracle because we 
need to ensure that the base of the prismoid, where x* is located, be a simplex. So, 
whenever a query hyperplane intersect the base, we m^e at most d ~ 2 additional 
queries to again make it a simplex. 

It be easily seen that this algorithm runs in linear time if the time complexity 
of oracle is 0(|5|). 


5.4 Intersection of Median of Straight Lines with 
a Given Line 


Let there be n non- vertical lines /j, where t =s in a Euclidean plane. Each 

of these lines is represented by a linear equation p = /,(i). The median level is 
defined as a function, F, such that F(a) is the ordinate of [n/2J-th intersection of 
the lines with the vertical line x = a. We can evaluate F(x) for any x in linear time 
by computing median of /,(x)’s. 

Let us compute intersection of median level and a given line /, which is represented 
by y = /(x). We are guarunteed a solution if the slope of this line is not the median 
slope if n is odd; does not lie between the [n/2j-th and [n/2 + Ij-th slope if n is 
even. 

For the sake of simplicity we assume that all the lines have positive slope and the 
line / is the x-axis. Then we have to determine x* such that F(x*) = /(«*) = y* = 0. 
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At X — + — oo and at x — y +oo function •F'(x) is same as the line of median slope. 
Therefore at x y +oo it is above the x-axis and at x — ► — oo it is below x-axis. 
Since function /^(x) is continuous, the equation F{x) = 0 has at least one solution. 
This problem was first posed by N. Megiddo to motivate the technique of paxametric 
searching. 

In the next subsections we briefly review how Megiddo and Cole applied their 
techniques to solve this problem. 


5.4.1 An Algorithm by Megiddo that Motivates the Tech- 
nique of Parametric Searching 

We first make an important observation. The order of lines changes only at inter- 
section points of a pair of lines. Hence, intersection points play an importamt role 
in the evaluation of x*. 

Let X,-,’ be the abscissae of the intersection point of lines /,• and Ij. 

We can compute x* naively as follows. First we identify intersection points of each 
pair of lines. Then we search for two values x^ and x^ such that F(x^) < 0 < 
F(x^) and there is no other Xy in the open interval (x^,x^). We search for such 
an interval by employing binaury search. This requires O(logn) F-evaluations and 
median computation of subsets of Xy whose cardinalities are N, N/2, JV/4, . . . points 
where N = Thus, this algorithm runs in O(ra^) time which is dominated by the 
evaluation of all the intersection points. 

Now we employ Megiddo’s technique to design an algorithm that runs in 0(n log^ n). 
We evaluate F(x*) with x* not specified by computing the median of /,’s. The 
median-finding algorithm compares values of function /,(x*) and fj{x*) for different 
z’s and ji’s. The outcome of sudi a comparison changes only at intersection point 
Xy . Hence, if F(xy ) > 0 then x* < Xy-; if F(xy ) < 0 then x* > Xy ; and if F(xy = 0 



100 


then a:* = Xij. Once we know the outcome of a comparison we proem! with the 
evaluation of median. Now the value of x* is restricted in the interval (~oo,x,ji) or 
(a;fj,+oo). At the end of median-finding algorithm we know the nredian line /, at 
X = x*. The intersection of this line with x-axis yields x\ Clearly, this algorithm 
is also 0(n®). Since the above algorithm is no better then the previous one we try 
another approach using AKS sorting networks. The AKS Sorting network employs 
P = 0(n) processors in every stage and sorts its input in O(logn) time. Instead of 
evaluating function F(x) at each of the Xj/s, we first wait for all the P values in 
a parallel iteration to be computed. Then we sort these P vduea and do a binary 
search to find the interval which contains x*. Then, we resolve ail the comparisons of 
this parallel iteration. "When all the iterations of AKS sorting network are completed 
then the lines axe sorted at x* and computation of x* is easy. 

Since each stage takes 0(n log n) time and there are O(logn) iterations, so the 
overall running time of this algorithm is 0(n log’ n). 


5.4.2 Algorithm Using Slowed Down Sorting Networks 

The above algorithm was further improved by Cole using slowed down sorting 
networks. 

A bri^ outline of his technique is as follows. A comparator in the sorting network 
is either active or inactive. An active comparator has its inputs determined but the 
outputs not determined. The other comparators axe inactive. An active comparator 
at depth d in the network is assigned a weight 1/4*^. The intersection point of padir 
of its input, Xij, is also assigned the same weight. 

In Megiddo’s algorithm we exhausted all the x./s of each iteration before going 
to the next iteration. Here we simultaneously deal will all the x<j’s of the active 
comparators with their respective weights. Accordingly, we evaluate F at weighted 
median of Xjj’s in the place of median, in the search of interval that contains 
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X*. Thus, we resolve comparisons in either the left half, (- 00 , x^], or the right 
half, {xto,oo) of the weighted median point These comparisons with weight 
approximately half the total weight become inactive and are replaced by active 
comparisons not more that one-fourth of the total weight as new weight. We reiterate 
these steps until we sort the lines at x*. We claim the following. 

Theorem 5.7 ([8]) At the completion ofkth iteration the active weight is bounded 
by (3A)*n/2, for Jb > 0. 

Theorem 5.8 ([8]) For k > b{d+ yjlogn), after kth iteration there are no active 
comparators at depth d, where d>0. 

Theorem 5.7 guarantees that in every iteration we drop a finite fraction of total 
weight whereas Theorem 5.8 guarantees that the algorithm terminates in O(logn) 
steps. Therefore, after C>(logn) steps we know the median of /i’s at x* and hence 
x“. The running time of this algorithm is 0(n log n). 

5.4.3 Application of Prune and Search Technique 

Now we apply prune and search technique to this problem. Let us refer to Algo- 
rithm 5 presented in the previous section. The input of this algorithm is set of lines 
in the plane. We do the following: 

• We prune the lines which are not the i:-th largest lines (initially k = [n/2j ) 
at X* and which lie above or below the trapezium given by Theorem 5.4. 

• When some of the lines are pruned we update value of k. 

• The oracle ft is constructed by evaluating F(j>) at x = p and determining 
position of x* with respect to p with the help of sign of F{p). 


This algorithm runs in lineaur time which is optimal. 
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5.5 Linear Programming 


The next problem that we solve is the linear programming problem in fixe<i dimen- 
sions. The linear programming problem is as follows. 

Minimise z 
Subject to 

(Ci, x) < z, 1 < t < n 


Esdh constraint defines a halfspace in a d-f 1 dimensional space. We try to compute 
the solution point x*. We again apply the Algorithm 5 for computing the fc-th 
largest hyperplane at x*. 

We construct the oracle 0, for this problem as follows. We recursively compute 
the optimum value constrained on the query hyperplane H. We can compute the 
relative position of the optimal point with respect to H using the gradient of largest 
constraint at the constrained optimal point on H. The base case in the recursion is 
a point P corresponding to 0-th dimension where the solution is trivial. 

We sort the constraints and choose the largest constraint at x*. Whenever a 
comparison is resolved the smaller hyperplane is pruned. Since the highest constraint 
is to be computed, Theorem 5.5 is applicable. This algorithm runs in linear time 
though it is not as efB.dent as that of Megiddo because of the large constants involved 
in the AKS sorting networks. 


5.6 Ham- sandwich Cuts 


Let Ai, where 1 < t < d, be d finite sets of points in The ham-sandwich cut of 
these is a hyperplane such that it bisects each of these sets. To compute this cut 
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we dualise the problem as usual. In the dual space the sets Ai transform to d sets 
of hyperplanes. We define this transformation as follows. 


Primal Space y = {m,x) + c {x,y) 

Dual Space Y = {x,X) Ay 


It should be noted that the above transformation conserves the incidence relation- 
ships. 


A ham-sandwich cut transforms to the point, say h, such that in eadi set there are 
exactly half the total number of hyperplanes that are above h. Consider the problem 
of computing h. If the median hyperplanes of each set at the point h are known 
then h can be computed by computing the intersection point of these. Hence, we do 
the following in the Algorithm 5. The elements that are to be sorted are the above 
hyperplanes and the ordering of the hyperplanes is to be done at h. A comparison 
in the sorting network would produce a “ridge”. We need to localise h with respect 
to these “ridges”. The comparisons corresponding to those ridges for which h is 
localised get resolved. 


In the oracle fl we need to know whidh side of a given hyperplane H does the point 
h lie. For this we use the method of Lo, Matousek and Steiger[17]. This provides 
us a method to determine if h lies inside a given infinite prism. This oracle Q first 
computes the median level of one of the sets on each vertical side of prism and 
then counts the total number of intersections below this level of the median levels of 
other d - 1 sets. If this number is odd then h lies inside the prism. The complexity 
of this algorithm is bounded by the complexity of constructing a fc*^-level in the 
d - 1 dimensional plane. So we get the same time bounds as in [17]. The linearly 
separated case in can also be solved similarly in linear time. 
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5.7 Centrepoints in d-dimensions 


In this section we discuss the problem of computing centrepoints of sets of points 
in d-dimensional space. Suppose 5 is a set of n points in The cenirr of S is 
union of all points p such that any closed halfspace contdning p contain.s at least 
fn/(d+ 1)] points. We obtjun a weaker type of centre, approximate or €-centre if 
we decrease the lower bound in the above definition of centre to fn(l — £)/(d-i- 1)] , 
where 0 < e < 1 . 

We first discuss the method of computing an approximate or e-centre of set S, where 

0 <e < 1. 

To compute an e-centrepoint we dualise the centrepoint problem. However, as 
opposed to ham-sandw'ich cut algorithm in the previous section, here we need to 
consider both the primal space as well as the dual space, in the algorithm. In 
particular, we search the centre in the dual space and prune the points in the primal 
space. As before we employ the technique by applying Algorithm 5. 

A centrepoint in the dual space will be a hyperplane (a centn-hyperplane) such 
that in every vertical line there would be at least fn/(d -f 1)] intersections of given 
hyperplanes above and below the intersection of centre-hyperplane. The dual of 
centre is union of such hyperplanes. The £-centre-hyperpIane is defined similarly. 

The elements to be sorted are the dual hyperplanes of points in S. We apply sorting 
on these hyperplanes. Since We do not execute the AKS sorting network to its 
completion, we can be vague about the sorting order and so we do not specify the 
exact point at which the sorting is taking place. We apply the algorithm just to 
partition the space such that each partition contains only small number of related 
unresolved comparisons. In the absence of localisation the oracle fl is not needed 
here. 

Whenever a hyperplane is queried in the Algorithm 5, we process the next steps 
allowing all the possible outcomes of the constrained search of the oracle In this 
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way whole of the region is partitioned into several sets each containing only small 
number of intersection with “ridges” related to unresolved comparisons. We iterate 
until the maximum of these numbers drops below Then we apply the 

prismoid method in each partition to compute 4£-centre. It can be seen that the 
region determined by fn/(d + 1)] lowest and highest levels is an 4£-centre. This 
is the consequence of Theorem 5.4. Clearly, each face of 4£-centre contains exactly 
fn/(d + 1)] points. 

Since we are not searching for a solution point, we do not prune points as usual. 
We provide another method of pruning below. In the earlier chapter we saw how 
pruning can be done in the plane for collection of three and four points. These 
points were respectively deleted and substituted by the Radon’s point. In fact, it 
ran be shown that (d+ l)r + 1 points in Si**, for any r, can be suitably substituted by 
their Tverberg’s point (a generalisation of Radon’s point) when certain conditions 
are met. Unfortunately, this also can not be applied in our case. To prune the 
points, we conjecture the following[33]. 

Conjecture 5.1 If there are (d + l)r + 9, where 0 < g < r points in such that 
each r-set is contained in some halfspace containing less than [n/(d+ 1)J points of 
S then these {d+l)r + q points can be substituted by q points. 

If this conjecture is proved then we will be able to discard at least a fraction of 
points in each iteration and a centrepoint in d dimensions can be found in linear 

time. 



Chapter 6 


Conclusions 


An optimal algorithm has been described in this thesis for computing the intersec- 
tion radius of a given set of planar objects containing points, straight lines, rays, 
line segments, wedges, half planes and planes. We have used a new approach of 
allowing addition of objects in the pruning process while ensuring a simultaneous 
net reduction of the input set. This provides one more example of the power and 
versatility of the prune and search paradigm. This algorithm has applications in 
many practical stabbing radii problems. 

We have also presented an optimal algorithm for computing a centrepoint for a 
planar finite set of points. This settles a long standing problem in computational 
geometry. It may be worthwhile to point out that an approximate centrepoint 
suffices in many applications but this exact algorithm, because of its efficiency, can 
very well replace the existing approximate centrepoint algorithms in two dimensions. 

A new technique has been presented to solve geometric optimisation problems using 
prune and search paradigm. We have been able to do this by the synthesis of paxa- 
metric searching and prune and search techniques. We have applied this tedbnique to 
compute ham-sandwich cuts and centrepoints in higher dimensions. This technique 
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has also been applied to the well known problem of linear programming. The results 
obtained match the best results extant for these problems. 


6.1 Further Resesirch Problems 


It would be worth investigating if there exist similar algorithms as presented in this 
thesis to compute intersection rauiius of objects in higher dimensions. The weighted 
version of this problem also remains unsolved. The intersection radius for a set 
of simple polygons cannot be computed using the same algorithm because of the 
inherent non-convexity of the distance function. 

Just as there are centrepoints, we can have other centre objects like centreline, 
centrecircle etc. The computation of these remain an open problem. The weighted 
version of the centrepoint problem also remains open. The other problems for a 
finite set of simple polygons in the plane to look at are area-centrepoints, perimeter- 
centrepoints and the like. The weighted centrepoint can be addressed as a special 
case of these problems. 

In this thesis, a centrei>oint of an explicit set of points is computed in linesur time. 
However, the centrepoint can not be computed within the same time bounds if the 
point set is implicit, such as the set of intersection points of straight line defined by 
any two points of a set of points in with a fixed plane. 

In the technique of parametric prune and search the open area of research is to 
provide efficient algorithms for pruning. The pruning technique described in this 
thesis has a doubly exponent factor of d in its time complexity which renders 
the algorithm impractical for higher dimensions. AKS sorting networks have also 
large constants of proportionality in their complexities. In practice, to solve the 
problem discussed in this thesis faster, randomised algorithms can be used. The 
other important question is how this technique can be applied to other parametric 
searching algorithms as well. It seems that the method presented in the thesis is a 
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particular case of Mataousek’s result[19]. It remains to be seen if this technique cam 
be applied to non-linear programming problems as well. 

F\irthermore, AKS sorting networks have built-in comparators that take two inputs 
each. The generalisation of comparison for points in higher dimensional space is their 
rotational order. So, if a network is designed which has such rotational ordering 
elements in plaice of comparators then it can be applied to many more problems 
using a similar technique. 
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